AD-AUO  016  NAVAL  POSTORAOUATE  SCHOOL  MONTCRCY  CA  F/9  4/2 

A  STATISTICAL  ANALYSIS  Of  MONTHLY  RAINFALL  FOR  MONTEREY  PENINSU— ETC I U> 
MAR  01  D  f  DAVIS 


Ifc  OLE  COPT  ADA110816 


LEW 

NAVAL  P08T6RADUATE  SCHOOL 

Monterey,  California 


A  STATISTICAL  ANALYSIS  OF  MONTHLY  RAINFALL  FOR 
MONTEREY  PENINSULA  AND  THE  CARMEL  VALLEY  IN 
CENTRAL  CALIFORNIA 

by 

David  Frederick  Davis 
March  1981 

Thesis  Advisor:  P.A.  Jacobs 

Approved  for  public  release;  distribution  unlimited 


i 


§2  02  11  C6S 


IICtfWTV  CLMUriCATlQH  QW  THlt  Mot  (W*m  Daw  toWfg 


REPODT  DOCUMENTATION  PAGE 


mad  msTitucrroNS 
BEFORE  COMPLETING  FORM 


NT'S  CATALOG  NUMICA 


4.  TITLE  CAR 4MMMI  I  TYPE  OF  REPORT  4  RERlOO  COVERED 

A  Statistical  Analysis  of  Monthly  Rainfall  Master's  Thesis; 

for  Monterey  Peninsula  and  the  Carmel  March  1981 _ 

Valley  in  Central  California  ••  rerformimo  oro.  report  numgsr 


David  Frederick  Davis 


».  performing  organization  name  ano  aooreii 

Naval  Postgraduate  School 
Monterey,  California  93940 


ii  controlling  office  name  ano  aoorem 

Naval  Postgraduate  School 
Monterey,  California 


RESVff  Affm>i 


NT  ft  m !• 


I  CMMJifAl  ottif) 


II.  HCHOHT  OATS 

March  1981 


II.  MUNICH  OR  HAOCI 

207 


II.  ICC  U  H|  TV  CL  All.  a  (AIM 

Unclassified 


l«.  OISTHIOuTION  STA 


Approved  for  public  release;  distribution  unlimited 


17.  OIITHICUTION  STATCNCNT  Cl  <*•  «D •trmt  m i—4  tn  f/Mt  20,  it  mttrnrmt  tmm 


I  fCNHw  m  mvmn  •*#»  H  on*  ttmrttr  Dr  D— «D 


Monthly  Rainfall  Models 
California  Rainfall 
Seasonal  Rainfall  Forecasting 
ARMA  Models 


Logistic  Analysis 
Logistic  Time  Series 


tfc.  aJMTHACT  (Cmttmm  m  mmtm  If  NMMHf  Nil  /HmUIT  Dr  D*#«D  MMN 

"^This  thesis  presents  a  statistical  analysis  of  the  monthly 
rainfall  for  the  Monterey  Peninsula  and  the  Carmel  Valley  in 
Central  California.  The  analysis  begins  with  the  simple  first- 
order  autoregressive  Markov  model,  which  is  found  to  be  weak. 
Next,  2x2  contingency  tables  are  used  to  identify  predictors, 
one  of  which  is  found  to  be  January  rainf all. ^Finally , 
logistic  analysis  is  used  to  quantify  the  prek^ctive  ability 


.  urng  in  a'I'tt  i  xn 


of  January. 

This  paper  attempts  to  analyze  rainfall  time  series  in  the 
statistical  sense.  No  attempt  is  made  to  provide  a  physical 
explanation  of  the  findings  from  the  point  of  view  of  a 
meteorologist . 


Accession  For 

"ntis  cnAfci 
DTIC  m 

U;*enno\n’t-d 

:  n. 


COPY 

.INSPECTED  . 


Bt - - - 

Dint  rib;;*.  ‘>n/ 

Avf.  U  !  •lit  .  Col  os 
Av awl /or 
Dint  D^c-ciol 


Approved  for  public  release;  distribution  unlimited 


A  Statistical  Analysis  of  Monthly  Rainfall  for 
Monterey  Peninsula  and  the  Carmel  Valley  in 
Central  California 


by 


David  Frederick  Davis 
Captain,  United  States  Army 
B.S.,  Colorado  School  of  Mines,  1972 


Submitted  in  partial 
requirements  for 


fulfillment  of  the 
the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 


)yO  /U  _ 

Deem  of  Information  and  Policy  Sciences 


3 


ABSTRACT 


This  thesis  presents  a  statistical  analysis  of  the 
monthly  rainfall  for  the  Monterey  Peninsula  and  the  Carmel 
Valley  in  Central  California.  The  analysis  begins  with  the 
simple  first-order  autoregressive  Markov  model,  which  is 
found  to  be  weak.  Next,  2x2  contingency  tables  are  used 
to  identify  predictors,  one  of  which  is  found  to  be  January 
rainfall.  Finally,  logistic  analysis  is  used  to  quantify 
the  predictive  ability  of  January. 

This  paper  attempts  to  analyze  rainfall  time  series 
in  the  statistical  sense.  No  attempt  is  made  to  provide 
a  physical  explanation  of  the  findings  from  the  point  of 
view  of  a  meteorologist. 
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I.  INTRODUCTION 


A.  THE  PROBLEM 

The  Monterey  F^.-r'Sula  Water  Management  District,  in 
Central  Califoi'iik  f  Astal  area  has  as  one  of  its  responsi¬ 
bilities  the  dur*.^  r-ti  recommend  and/or  impose  water  rationing 
on  its  constituents.  To  do  this  in  a  rational  way  requires 
the  District  to  have  some  formula  for  predicting  future  water 
availability.  Although  the  techniques  of  modern  meteorology 
are  becoming  more  sophisticated  and  exact  there  is  still  the 
inability  to  make  good  long-range  predictions.  This  thesis 
analyzes  three  series  of  Monterey  County  monthly  rainfall 
data  by  purely  statistical  methodology  in  order  to  identify 
possible  predictive  formulas. 

B .  NOTATION 

Rainfall  will  be  denoted  as  Rfc  m  which  will  represent 
inches  of  rain  recorded  for  the  t^  year  and  the  m^*1  month. 
The  year  to  be  used  is  the  California  Water  Year  which  begins 
in  October  and  ends  the  following  September.  Thus  R^  ^  is 
the  monthly  rainfall  for  October  of  year  ’ 1 1  and  Rg  g  is 
the  monthly  rainfall  of  May  of  year  '6'. 

An  overstruck  bar  as  in  R. .  will  indicate  the  arithmetic 
average  of  a  variable;  in  this  case  it  is  the  arithmetic 
average  of  all  years  and  months  of  rainfall.  is  the 

average  of  rainfall  over  the  years  for  month  m;  Rt# 


represents  the  yearly  average  for  year  t. 


C.  METHODS  OF  ANALYSIS 

Three  methods  were  used  to  analyze  the  data.  The  first 
method  was  to  model  the  series  using  autoregressive  moving 
averages  as  described  in  Box  and  Jenkins  [Ref.  1],  The 
second  was  to  use  2x2  contingency  tables  to  identify  possible 
predictors.  The  third  was  logistic  regression  to  quantify 
the  findings  of  the  2x2  contingency  table  analysis.  These 
three  methods  will  be  described  in  further  sections  of  this 
paper. 

1.  ARMA(p,q)  Models 

A  widely  used  approach  to  time  series  modeling  pro¬ 
posed  by  Box  and  Jenkins  is  the  ARMA(p,q)  model.  This  model 
is  actually  a  joining  of  two  types  of  model,  the  autoregres¬ 
sive  and  the  moving  average. 

In  the  noation  of  Box  and  Jenkins: 
let  {Zfc,  t»l,2,...,n}  be  a  time  series,  then  am  ARMA (p ,q) 
process  may  be  written  as: 

X  f  X 

Zt  »  +  ...  +  0pZt_p  +  at  -  9 1at_1  -  •••"dqat-q 

the  {afc  ,  t»l,2,...,n}  are  assumed  to  be  random  shocks 
distributed  as  independent  and  identically  distributed  (iid) 
random  variables  with  meam  zero  and  variance  a*  amd 

a 

Z  «  Z.  -  Z. .  The  further  assumption  of  normality  is  also 

w  w 

usually  made. 

For  purposes  of  this  paper,  a  mapping  of  Rt  m 
,  r«12(t-l)+m  was  made,  and  an  ARMA  analysis  was 


into  Z 


then,  conducted  on  this  index  transformed  series.  This 


analysis  is  described  in  section  III. 

2 .  2x2  Table  Analysis 

In  the  validation  of  section  IV  it  is  found  that 
the  ARMA  model  is  not  very  successful  in  describing  the 
data.  In  section  V  the  data  is  analyzed  by  means  of  2x2 
contingency  tables.  These  tables  are  good  tools  for  explo¬ 
ratory  data  analysis  in  that  they  provide  a  visual  display 
of  the  data.  Statistical  procedures  based  on  the  null 
hypothesis  of  independence  can  be  used  to  quantify  the 
departure  from  independence.  The  theory  of  2x2  tables, 
and  contingency  tables  in  general  may  be  found  in  Fleiss 
[Ref.  3],  Dixon  and  Massey  [Ref.  5],  Brownlee  [Ref.  6], 
and  Mood,  Graybill,  and  Boes  [Ref.  7]. 

For  this  paper,  the  contingency  table  approach  is 
used  to  identify  a  month  or  group  of  months  of  a  year  whc.a 
rainfall  cam  serve  as  a  predictor  for  the  rainfall  during 
the  remaining  months  of  the  year.  One  predictor  that  was 
suggested  is  the  rainfall  in  the  month  of  January. 

3 .  Logistic  Analysis 

Once  a  predictor  is  tentatively  identified  it  becomes 
necessary  to  quantify  the  degree,  direction  and  accuracy  of 
the  predictor. 

A  logistic  amalysis  is  conducted  by  dividing  the  data 
for  a  year  into  two  3ets,  the  predictor  or  control  set,  and 
predictand  or  complement  set.  For  this  analysis,  the  predictor 
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is  the  logged  anomaly  of  January  rainfall  for  the  year; 


that  is,  if  denotes  the  predictor  or  control  for  year 
t ,  then 


N 


Xt  >  £n(Rt(4) 


1  ^ 

-  -  y 

"  t»l 


£n(Rt/4} 


1.2 


(The  logarithm  is  used  to  better  symmetrize  the  model.)  The 
complement  is  the  raw  anomaly  of  the  total  rainfall  for  the 

I 

immediately  subsequent  eleven  months;  that  is,  if  Y  denotes 
the  complement  for  year  t,  then; 


12 

(  l  Rt 
m=5  ' 


m  +  JlRt+1'm) 


1.3 


.  N-l  12  3 

*  H=T  I.  <  l  Rt,m  +  Vl,m> 

t=l  m=5  '  m=l 

Finally,  the  data  are  transformed  into  a  binary  representa¬ 
tion,  relative  to  zero  as; 


1.4 


0 

if  Y 1 1  < 

0 

1 

if  Y*  > 

0 

In  section  VI  the  model  fit  is 


P ( Y=1 i X=x)  = 


.0+3  x 


1  +  e 


a+3 


1.5 


Where  x  is  as  before  and  P(Y=l|X=x)  is  interpreted  as: 

"the  conditional  probability  that  the  subsequent  eleven  month 
total  rainfall  will  be  above  its  mean,  given  that  the  logged 
anomaly  of  January  rainfall  was  ’x’". 
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II.  THE  DATA 


i 


A.  GENERAL 

Three  data  sets  were  used  for  this  analyis.  The  location 
at  which  these  data  sets  were  gathered  is  shown  in  Figure  1. 

As  the  figure  indicates,  two  of  the  data  sets  are  on  the 
Monterey  Peninsula  proper,  while  the  third  set,  SC,  represents 
the  Carmel  River  Watershed  at  the  San  Clemente  Dam. 

Although  data  exists  in  all  cases  to  the  present,  all 
three  sets  were  truncated  at  September  of  1974.  The  remain¬ 
ing  data,  up  through  September  of  1980  was  reserved  for 
validation  of  the  models  and  methodology. 


The  data  coordinates  are : 


Data 

set 

RN: 

36o 

35* 

42" 

North  Latitude 

121° 

54' 

43" 

West  Longitude 

Data 

set 

FL: 

36° 

35  • 

30" 

North  Latitude 

121° 

56  ' 

30" 

West  Longitude 

Data 

set 

SC: 

36° 

26' 

12” 

North  Latitude 

121° 

42  • 

30" 

West  Longitude 

Monterey  Bay 


Figure  1.  Location  of  rainfall  data  sets 
and  the  years  available. 


B.  DATA  SET  RN 


Data  set  RN  consists  of  monthly  rainfall  amounts  gathered 
by  Professor  R.J.  Renard,  Cooperative  Observer  for  the  National 
Weather  Service  Climatological  Station,  Monterey,  California. 
The  data  set  begins  in  June  1951  and  currently  terminates  in 
September  1980.  As  was  stated  above,  the  analysis  was  con¬ 
ducted  only  on  that  data  between  and  including  October  1951 
and  September  1974. 

1.  Raw  Data 

Appendix  A  contains  a  listing  of  data  set  RN.  Figure 
2  shows  the  raw  data  set.  Month  1  is  October  1951,  month  148 
is  January  1964,  and  up  to  month  288  which  is  Seti^Aiber  j.^',4. 

As  cam  be  seen  the  data  are  strongly  seasonal.  This  is  enough 
to  indicate  that  the  series,  as  stated,  is  highly  non¬ 
stationary. 

The  data  presented  so  far  deals  with  only  monthly 
data.  Next  to  be  considered  is  the  series  of  yearly  total 
rainfalls.  The  results  are  shown  in  Figure  3  (Yearly  total 
rainfall) ,  4  (Correlogram  of  yearly  rainfall) ,  and  Table  1 
(Estimated  Autocorrelations).  In  this  case,  the  correlogram 
indicates  stationarity  and  independence  of  the  yearly  series. 

A  plot  of  the  lag  one  relationships,  Figure  5,  reinforces 
this  indication  of  independence. 
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0.000 


Figure  2.  Monthly  rainfall  in  inches  for  data  set  RN 


The  ccrrelograms  and  Partial  Correlograms  to  follow 
indicate  the  95%  approximate  significance  levels  using 
dashed  lines.  For  development  of  these  significance 
levels  see  Box  and  Jenkins  [Ref.  1]. 


Figure  3.  Yearly  total  rainfall  for  data  set  RN 
(1951  -  1974)  . 
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Figure  4.  Correlogram  of  Yearly  total  rainfall 
for  data  set  RN. 


TABLE  1 

ESTIMATED  AUTOCORRELATIONS  OF  YEARLY 
TOTAL  RAINFALL  FOR  DATA  SET  RN 


AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

-.200 

7 

.044 

2 

.  109 

8 

-.109 

3 

-.295 

9 

-.009 

4 

.186 

10 

-.139 

5 

-.249 

11 

.248 

12 

-.211 

29 
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monthly  data  points 
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II. 6 


These  statistics  for  data  set  RN  are  shown  in  Table  2,  and 
illustrated  in  Figure  6.  In  the  same  way  as  the  raw  data 
mapped  into  a  series,  a  series  is  created  from 


Zr  »  Rfc  m  ,r=12 (t-1) +m 


II. 7 


TABLE  2 

MONTHLY  MEANS  AND  VARIANCE  FOR  DATA  SET  RN 


MONTH 

MEAN 

VARIANCE 

1 

.677 

.5292 

2 

2.478 

3.6354 

3 

3.355 

5.9147 

4 

4.124 

5.4146 

5 

2.592 

4.5122 

6 

2.6  39 

3.7414 

7 

1.  708 

2.7699 

8 

.434 

.2494 

9 

.219 

.1004 

10 

.058 

.0056 

11 

.104 

.0129 

12 

.275 

.  3953 
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Figure  7.  Monthly  rainfall  anomalies  in  inches 
for  data  set  ?N 
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Figure  8.  Correlogram  of  the  monthly  rainfall 
anomalies  for  data  set  RN 


TABLE  3 

ESTIMATED  AUTOCORRELATIONS  OF  MONTHLY  RAINFALL 
ANOMALIES  FOR  DATA  SET  RN 


LAG 

VALUE 

LAG 

VALUE 

1 

.249 

14 

-.057 

2 

-.090 

15 

-.109 

3 

-.059 

16 

-.006 

4 

.041 

17 

.  067 

5 

.032 

18 

.004 

6 

-.035 

19 

.013 

7 

-.011 

20 

-.007 

f 

8 

-.043 

21 

.063 

9 

-.073 

22 

.023 

10 

-.091 

23 

-.  066 

11 

.076 

24 

-.044 

12 

-.020 

25 

.021 

13 

-.012 
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3.  Logged  and  Swept  Data 

The  data  should  now  be  stationary  in  the  means. 
However,  as  seen  in  Table  2,  the  variances  of  monthly 
rainfall  amounts  are  not  homogeneous.  Kilmartin  [Ref.  10] 
discusses  various  transformations  of  the  data  to  remove 
this  heteroskedacity .  A  plot  of  the  variance  versus  mean. 
Figure  9  below,  indicates  that  the  logarithmic  transform 
of  the  data  might  be  useful. 
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Figure  9.  Plot  of  monthly  variance  against 
monthly  means  for  data  set  RN 
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Since  the  data  contain  zeros,  the  following  modified 
logarithmic  transformation  is  done 


II. 8 


where  the  effect  of  the  addition  of  the  one  is  mostly  to 

preserve  the  mapping  of  zeros  into  zeros.  A  more  in  depth 

discussion  of  this  transformation  is  found  in  Kilmartin. 

The  mapping  is  performed  again  as  before  and  R'  and 

.  m 

S'  are  calculated  in  a  manner  similar  to  II. 6  and 

.  m 

shown  in  Table  4  and  Figures  10  and  11. 


TABLE  4 

MONTHLY  MEANS  AND  VARIANCE  FOR 
LOGGED  DATA  SET  RN 


MONTH 

MEAN 

VARIANCE 

1 

.438 

.1549 

2 

1.092 

.3454 

3 

1.312 

.  3538 

4 

1.539 

.2003 

5 

1.104 

.3795 

6 

1. 133 

.  3706 

7 

.  854 

.2728 

8 

.319 

.0767 

9 

.176 

.0401 

10 

.054 

.0045 

11 

.094 

.0094 

12 

.185 

.0849 
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These  transformations,  the  logarithm  followed  by 
the  removal  of  the  monthly  means  of  the  logged  data,  result 
in  the  series  listed  in  Appendix  A  and  described  in  Figures 
12  and  13  with  Table  5. 

These  displays  indicate  that  a  suitably  stationary 
series  has  been  obtained.  Other  methods,  such  as  differencing, 
scaling,  and  Box-Cox  transformations,  see  Hipel  [Ref.  11], 
were  tried  but  with  less  success. 
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Figure  12a.  Logged  anomalies  of  monthly 
rainfall  for  data  set  RN.  Months  1  -148 
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Figure  12b.  Logged  anomalies  of  monthly 
rainfall  for  data  set  RN 
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Figure  13.  Correlogram  of  logged  anomalies  of 
monthly  rainfall  from  data  set  RN 
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TABLE  5 


ESTIMATED  AUTOCORRELATIONS  OF  LOGGED 
ANOMALIES  OF  MONTHLY  RAINFALL  FROM 
DATA  SET  RN 


AUTOCO RRE LAT I ON S 


LAG 

VALUE 

LAG 

VALUE 

1 

.191 

14 

-.033 

2 

-.102 

15 

-.084 

3 

-.095 

16 

-.024 

4 

.071 

17 

.062 

5 

.056 

18 

-.015 

6 

-.053 

19 

-.008 

7 

-.009 

20 

.013 

8 

-.014 

21 

.013 

9 

-.022 

22 

-.038 

10 

-.069 

23 

-.024 

11 

.024 

24 

-.012 

12 

-.004 

25 

.032 

13 

.032 

C.  DATA  SET  FL 

The  label  for  these  data  derives  from  its  location, 
Forest  Lake,  on  the  Monterey  Peninsula,  in  Pebble  Beach, 
California.  Data  set  FL  consists  of  monthly  rainfall 
figures  gathered  by  the  California-American  Water  Company 
since  1896.  Although  this  data  set  started  quite  early, 
the  data  prior  to  1937  has  frequent  missing  observations. 
Therefore,  this  data  set  is  taken  as  October  1937  through 
September  1974,  with  October  1974  through  September  1980 
reserved  for  validation. 

Analysis  of  this  data  set  is  identical  to  that  of  data 
set  RN,  therefore  only  the  pertinent  figures  and  tables 
are  shown . 
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Figure  14a. 


Months, 1  -  296  of  rainfall  in 
inches  for  data  set  FL 
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Figure  14b.  Months  297-444  of  rainfall 
in  inches  for  data  set  FL 
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Figure  15  .  Yearly  total  rainfall  for  data  set  FL 
(1937  -  1974) . 
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Figure  16.  Correlogram  of  yearly  total  rainfall 
for  data  set  FL 


TABLE  6 

ESTIMATED  AUTOCORRELATIONS  OF  YEARLY  TOTAL  RAINFALL 

FOR  DATA  SET  FL 

AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.010 

14 

-.028 

2 

.099 

15 

.105 

3 

-.022 

16 

-.087 

4 

.207 

17 

.085 

5 

-.214 

18 

-.018 

6 

.028 

19 

-.245 

7 

-.003 

20 

-.237 

8 

-.061 

21 

-.035 

9 

-.138 

22 

-.031 

10 

-.060 

23 

-.106 

11 

.236 

24 

-.096 

12 

-.119 

25 

.120 

13 

-.151 
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TABLE  7 


Swept  Data 


MONTHLY  MEANS  AND  VARIANCE  FOR  DATA  SET  FL 


MONTH 

MEAN 

VARIANCE 

1 

.744 

.  4895 

2 

2.235 

3.7444 

3 

3.049 

4.2480 

4 

3.537 

4.4240 

5 

2.999 

5.9492 

6 

2.724 

3.1743 

7 

1.559 

2.4505 

8 

.449 

.2033 

9 

.153 

.0417 

10 

.077 

.0081 

11 

.115 

.0081 

12 

.189 

.1350 
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Figure  17.  Monthly  means  for  data  set  FL 


44 


RRINfftLL  HNOMflLILS  IN  INCHES  'tJ  *  RHINE RLL  RNOMRLICS  IN  INCHLS 


DftTn  IS  .,i£^F‘; 


MONTHS  OF  ZRTR 


Figure  18a.  Months  1  -  296  of  rainfall  anomalies 
in  inches  for  data  set  FL 
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Figure  18b.  Months  297-444  of  rainfall  anomalies 
in  inches  for  data  set  FL 
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Figure  19.  Correlogram  of  monthly  rainfall 
anomalies  for  data  set  FL 
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TABLE  8 


ESTIMATED  AUTOCORRELATIONS  OF  MONTHLY 
RAINFALL  ANOMALIES  FOR  DATA  SET  FL 


AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.244 

14 

-.053 

2 

-.007 

15 

-.043 

3 

-.056 

16 

-.002 

4 

.027 

17 

.039 

5 

.016 

18 

-.014 

6 

-.026 

19 

-.011 

7 

-.022 

20 

-.031 

8 

-.021 

21 

.023 

9 

-.051 

22 

.067 

10 

-.020 

23 

-.007 

11 

.077 

24 

.041 

12 

.059 

25 

.039 

13 

.068 
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Figure  20.  Plot  of  monthly  variance  against  monthly 
means  for  data  set  FL 
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TABLE  9 


MONTHLY  MEANS  AND  VARIANCE  FOR  LOGGED  DATA  SET  FL 


MONTH 

MEAN 

VARIANCE 

1 

.484 

.  1440 

2 

1.007 

.  3440 

3 

1.269 

.2781 

4 

1.398 

.2579 

5 

1.218 

.  3421 

6 

1.184 

.  3006 

7 

.796 

.2717 

8 

.  338 

.0589 

9 

.  130 

.0234 

10 

.071 

.0061 

11 

.106 

.0064 

12 

.146 

.0044 
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L06GCD  *'Er>NS 


Figure  21.  Monthly  means  of  logged 
data  set  FL 
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Figure  22. 


Plot  of  monthly  variance  against  monthly 
means  for  logged  data  set  FL 
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Figure  23a.  Months  1  -  148  of  logged  rainfall 
anomalies  for  data  set  FL 
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Figure  23b.  Months  149  -  444  of  logged  rainfall 
anomalies  from  data  set  FL 
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Figure  24.  Correlogram  of  logged  anomalies  of  monthly 
rainfall  from  data  set  FL. 


TABLE  10 


ESTIMATED  AUTOCORRELATIONS  OF  LOGGED 
ANOMALIES  FROM  MONTHLY  RAINFALL  OF 
DATA  SET  FL 


AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.185 

14 

-.052 

2 

-.020 

15 

-.021 

3 

-.081 

16 

-.020 

4 

.046 

17 

.040 

5 

.043 

18 

-.040 

6 

-.050 

19 

-.025 

7 

-.024 

20 

-.014 

3 

-.010 

21 

.007 

9 

-.027 

22 

.019 

10 

.004 

23 

-.015 

11 

.031 

24 

.076 

12 

.084 

25 

.047 

13 

.074 
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D.  DATA  SET  SC 


The  label  for  this  data  derives  for  its  location,  San 
Clemente  Dam,  on  the  Carmel  River  in  Central  Calif ronia, 
approximately  26  kilometers  southeast  of  data  sets  RN  and 
FL  on  the  Monterey  Peninsula.  Data  set  SC  consists  of 
monthly  rainfall  figures  gathered  by  the  California-American 
Water  Company  since  1926. 

Analysis  of  this  data  set  is  again  very  close  to  that 
of  the  previous  data  sets  and  only  the  displays  will  be 
given. 

1.  Raw  Data 
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Figure  25a.  Months  1  (October  1926)  -  148 
(January  1938)  of  rainfall  in  inches  for 
data  set  SC. 
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Fiugre  25b.  Months  149  -  444  of  rainfall 
for  data  set  SC 
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Figure  25c.  Months  445  -  576  of  rainfall  for 
data  set  SC 
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Figure  26.  Yearly  total  rainfall  for  data 
set  SC  (1926  -  1974) 
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Figure  27.  Correlogram  of  yearly  total 
rainfall  for  data  set  SC 


TABLE  11 

ESTIMATED  AUTOCORRELATIONS  OF  YEARLY 
TOTAL  RAINFALL  FOR  DATA  SET  SC 

AUTOCORRELATIONS 


LAG 

VALUE 

1 

-.050 

2 

.  135 

3 

-.158 

4 

.16  8 

5 

-.042 

6 

.081 

7 

-.084 

8 

-.217 

9 

-.111 

10 

-.107 

11 

.158 

12 

-.260 

13 

-.110 

LAG 

VALUE 

14 

-.109 

15 

.161 

16 

.077 

17 

.116 

18 

.025 

19 

-.214 

20 

.066 

21 

.019 

22 

-.028 

23 

-.101 

24 

-.030 

25 

-.042 
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TABLE  12 


MONTHLY  MEANS  AND  VARIANCES 
FOR  DATA  SET  SC 


MONTH 

MEAN 

VARIANCE 

1 

.698 

.5945 

2 

2.175 

4.2382 

3 

3.940 

10.1783 

4 

4.599 

8.4899 

5 

4.353 

13. 3443 

6 

3.080 

5.2744 

7 

1.700 

3.4486 

8 

.431 

.1912 

9 

•  111 

.0451 

10 

.017 

.0055 

11 

.037 

.0125 

12 

.103 

.1040 
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Figure  28.  Monthly  means  for  data  set  SC 


13  FFor? 


1 1 . 000  , 


* 


z 


X 

x 

o 

z 

X 


X 


*  . 

\  '  \  .  ,  ..  ; 

»«*.  . . |i|wmipi4A(i«mH« 


X 

& 


i  ' 


V  I  »>  »  - 

I 


I 

1/ 


-4.343 

443 


MONTHS  CF  DHTfl 


592 


Figure  29b.  Months  297  -  576  anomalies  in 
inches  for  data  set  SC 


HUTOCORRELHT IONS 


CORRELOGFflM 


1 .0  r 


a.5  h 


0.0 


-0.5 


V 


-1 .0  L 


25 


LAG 


Figure  30.  Correlogram  of  monthly  rainfall 
anomalies  for  data  set  SC 


TABLE  13 

ESTIMATED  AUTOCORRELATIONS  OF  MONTHLY 
RAINFALL  ANOMALIES  FOR  DATA  SET  SC 

AUTOCORRELATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.140 

14 

-.088 

2 

-.039 

15 

-.051 

3 

-.021 

16 

-.006 

4 

.012 

17 

.015 

5 

-.001 

18 

-.006 

6 

-.019 

19 

-.008 

7 

.003 

20 

-.027 

8 

-.013 

21 

-.003 

9 

-.038 

22 

-.122 

10 

-.038 

23 

-.006 

11 

.014 

24 

.011 

12 

.051 

25 

.011 

13 

.102 
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3 .  Logged  and  Swept  Data 


ME  fir-13  V= 

M.000  r 


,  fir:  I  fif C-  Z"Z 


K 


K 


0. 000 


0.000 


MCRNS 


1 4 . 000 


Figure  31.  Plot  of  monthly  variance  against  anomally 
means  for  data  set  SC. 


TABLE  14 

MONTHLY  MEANS  AND  VARIANCES  OF 
LOGGED  DATA  SET  SC 


MONTH 

MEAN 

VARIANCE 

1 

.444 

.1623 

2 

.961 

.  3999 

3 

1.408 

.  3949 

4 

1.583 

.3117 

5 

1.444 

.4928 

6 

1.243 

.  3556 

7 

.817 

.  3247 

8 

.  320 

.0726 

9 

.092 

.0227 

10 

.015 

.0040 

11 

.031 

.0083 

12 

.075 

.0351 
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Figure  32.  Monthly  means  of  logged  data  set  SC 
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Figure  33.  Plot  of  monthly  variance  against  monthly 
means  for  logged  data  set  SC 
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Figure  34a. 


Months  1  -  296  of  logged  rainfall 
anomalies  from  data  set  SC 
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Figure  34b.  Months  297  -  576  of  logged  rainfall 
anomalies  from  data  set  SC 
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Figure  35.  Correlogram  of  logged  anomalies  of 
monthly  rainfall  from  data  set  SC 


TABLE  15 

ESTIMATED  AUTOCORRELATION  OF  LOGGED 
ANOMALIES  OF  MONTHLY  RAINFALL  FROM 
DATA  SET  SC 

AUTOCORRE  LATIONS 


LAG 

VALUE 

LAG 

VALUE 

1 

.096 

14 

-.057 

2 

-.065 

15 

-.022 

3 

-.066 

16 

-.023 

4 

.038 

17 

.021 

5 

.012 

18 

-.017 

6 

-.061 

19 

-.005 

7 

.023 

20 

-.008 

8 

-.001 

21 

-.011 

9 

-.056 

22 

.092 

10 

-.021 

23 

-.019 

11 

.007 

24 

.050 

12 

.091 

25 

.002 

13 

.091 
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III.  FIRST  ORDER  MARKOV  MODEL 


A.  THEORY 

As  first  shown  by  equation  1.1,  the  general  ARMA (p , q) 
model  is: 


zt  *  *1  zt-i  *■ 


•+*p  Zt-P  -  ei  at-i 


-  .  .  .  -0. 


t-q 


III.  1 


The  development  and  discussion  of  this  type  of  model  is 
contained  in  detail  in  Box  and  Jenkin  [Ref.  1]  and  Nelson 
[Ref.  8].  The  modeling  process  is  a  three  fold  procedure. 
The  parts  are: 

(1)  Identification 

(2)  Estimation 

(3)  Diagnosis. 

Identification  is  conducted  using  the  correlogram  and 
a  plot  of  the  partial-autocorrelations  (or  partial  correlo¬ 
gram)  .  The  partial  autocorrelations  are  related  to  the 
autocorrelations,  see  Box  and  Jenkins  [Ref.  1],  Nelson 
[Ref.  8],  or  Richards  and  Woodall  [Ref.  12 J .  These  partial 
autocorrelations  are  used  to  determine  the  order  of  the 
moving  average  process  much  like  the  autocorrelations  may 
be  used  to  determine  the  order  of  the  auto-regressive 
process . 

Once  the  autocorrelations  and  partial  autocorrelations 
have  been  found,  the  degree  of  the  ARMA  may  be  estimated  by 
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techniques  described  in  Box  and  Jenkins,  Nelson  or  Richards 
and  Woodall.  Each  of  the  data  sets,  once  logged  and  swept, 
indicated  that  the  most  probable  model  was  an  ARMA(1,0)  or 
AR ( 1)  or  more  commonly  a  first-order  autoregressive  Markov 
model.  This  model  is  simply; 

K  *  P^t-i  +  at  m.2 

where  the  p  is  the  autocorrelation  of  lag  one.  Thus,  this 
model  indicates  that  any  persistence  in  the  data  are  condi¬ 
tionally  independent  of  the  past  given  the  lag  one  value. 

Subsections  B,  C,  and  D  below  show  this  model  as  applied 
to  the  three  data  sets  of  interest.  The  residuals  of  the 
model  2fc  -  pZt_^  are  examined.  The  residuals  appear  to 
be  independent,  however,  they  do  not  appear  to  be  normally 
distributed;  for  example,  there  is  a  high  peak  around  zero. 

One  possible  reason  for  this  discrepancy  may  be  the  dichotomy 
of  winter  and  summer  rain  as  indicated  in  Tables  2,  4,  7,  9, 
12,  and  14.  The  existence  of  months  with  zero  rainfall  during 
the  summer  suggests  that  one  should  consider  the  summer,  when 
rain  is  sparse,  completely  separate  from  the  winter  when  rain 
is  more  abundant.  Therefore,  also  shown  in  the  subsections 
below  is  the  autoregressive  model  applied  to  the  winter  months 
only.  This  is  accomplished  by  stripping  out  months  9  through 
12  (June  through  September)  of  the  data  sets  and  treating  the 
remaining  data  as  a  continuous  set.  In  other  words,  the 
first  ten  months  are  then 

Rl,l'  R1 , 2 '  R1,3'R1,4'R1,5'  R1,6'R1,8'R2,1'R2,2‘ 
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The  appropriate  correlograras .and  partial  correlograms  are 
displayed  prior  to  the  model  applications. 

B.  DATA  SET  RN 

1 .  Twelve  Month  Series 

This  data  set  is  described  in  section  II. b.  The 
remaining  diagnostic  device  needed  is  the  partial  correlogram 
of  Figure  36  and  the  corresponding  values  in  Table  16. 
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Figure  36.  Partial  correlogram  of  the  logged 
rainfall  anomalies  of  data  set  RN 
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TABLE  16 


ESTIMATED  PARTIAL-AUTOCORRELATIONS  FOR 
LOGGED  RAINFALL  ANOMALIES  OF  DATA  SET  RN 


LAG 

VALUE 

LAG 

VALUE 

1 

.191 

14 

-.040 

2 

-.  144 

15 

-.073 

3 

-.047 

16 

.001 

4 

.092 

17 

.054 

5 

.006 

18 

-.067 

6 

-.058 

19 

.039 

7 

.037 

20 

.012 

8 

-.034 

21 

-.001 

9 

-.028 

22 

-.043 

10 

-.056 

23 

.013 

11 

.048 

24 

-.043 

12 

-.041 

25 

-.034 

13 

.047 

The  model  of  interest  is  then 

i't  *  .19125^  +  afc  III.  3 

where  the  random  shocks  {afc}  are  assumed  to  be  distributed 
iid  N(0,c2)  and  a2  is  estimated  as 

d  <3. 

■TT  X  <5t  *  •19lit-l)  111  ■* 

The  goodness  of  this  fit  may  be  viewed  in  two  ways. 
Firstly,  are  the  residuals,  {at>  independent?  Secondly, 
are  the  residuals  distributed  as  Normal  (Gaussian;  random 
variables?  A  plot  of  the  residuals  follows  in  Figure  37. 

The  question  of  independence  is  addressed  in  Figure  38 
(Correlogram) ,  Figure  39  (Lag  one  plot) ,  Figure  40  (Residuals 
vs.  lag  one) ,  and  Table  17  (Turning  points) .  For  a  discussion 
of  the  usefulness  of  the  turning  points  see  Kendall  [Ref.  14]. 
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All  of  these  displays  and  tests  tend  to  indicate  that  the 
residuals  are  in  fact  serially  independent.  The  statistics 
of  the  residuals  are  in  Table  18.  A  Normal  Plot  of  the 
residuals  (Figure  41)  ,  in  which  the  sample  is  normalized  by 
removing  the  mean  and  scaling  by  the  standard  deviation  and 
then  plotted  on  normal  paper,  should  yield  a  nearly  straight 
line  corresponding  to  the  dashed  line  of  the  figure.  The 
Normal  Plot  accompanied  by  the  sample  histogram  (Figure  42) 
addresses  the  normality  of  these  data.  As  may  be  seen  from 
the  kurtosis,  the  fluctuations  of  the  sample  CDF  near  the 
midpoint,  and  the  peak  of  the  histogram,  the  normality  of 
this  data  are  questionable.  To  confirm  this  a  chi-squared 
goodness  of  fit  test  was  conducted  yielding  a  value  of  49.18 
with  17  degrees  of  freedom,  again  rejecting  any  hypothesis 
of  normality  at  a  significance  level  of  5xl0~5. 
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Figure  37.  First  order  Markov  residuals  from  logged 
rainfall  anomalies  of  data  set  RN .  Months  149  -  292 
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Figure  38.  Auto  correlations  of  residuals  from 
first  order  Markov  process  applied  to  the  logged 
rainfall  anomalies  of  data  set  RN 


-1.141  1.246 

LAGGED  ELEMENTS 
W<  t - 1  ) 

Figure  39.  Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  RN 
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Figure  40.  First  order  Markov  residuals  versus 
lag  one  data  point  from  logged  rainfall  anomalies 
of  data  set  RN 


TABLE  17 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING  POINTS 
AND  ACTUAL  AND  EXPECTED  PHASE  FREQUENCIES 
FOR  THE  FIRST  ORDER  MARKOV  RESIDUALS  FROM 
THE  LOGGED  RAINFALL  ANOMALIES  OF  DATA  SET  RN 


NUMBER 

OF  TURNING  POINTS 

=  191 

E  [P  ]  * 

190.667  V[P] 

=  15.899 

PHASE 

LENGTHS 

D 

OBS  . 

E[*] 

1 

117 

118.8 

2 

56 

52.1 

3 

15 

14.9 

4 

3 

3.2 

5 

0 

.6 

6 

0 

.1 

7 

0 

0.0 
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TABLE  18 


GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 
ANOMALIES  OF  DATA  SET  RN 


Moments 


Mean 

-.001 

Variance 

.117 

Skewness 

-.066 

Kurtosis 

.523 

Percentiles 

Minimum 

-1.141 

Lower  Sixteenth 

-  .745 

Lower  Eight 

-  .463 

Lower  Quartile 

-  .174 

Median 

-  .014 

Upper  Quartile 

.211 

Upper  Eight 

.4  36 

Upper  Sixteenth 

.706 

Maximum 

1.246 

Figure  41.  Standardized  normal  plot  of  first 
order  Markov  residuals  from  logged  rainfall 
anamlies  of  data  set  RN 
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HISTOGRAM  :r 

riRST  ORDER  MfiKOV  RESIDUALS  ' 

Figure  42.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  RN 

2 .  Winter  Series 

As  stated  above,  the  number  of  summer  months  with 
zeros  indicated  that  a  look  at  the  winter  months  only  might 
be  worthwhile.  The  Figures  43  (Winter  months) ,  44  (Correlo- 
gram) ,  and  45  (Partial  correlogram) ,  which  deal  only  with 
winter  months,  still  indicate  a  first  order  autoregressive 
Markov  model  as; 

*  .21825^  +  at-;L  III. 5 
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Figure  43.  Winter  months  only  of  logged  rainfall 
anomalies  of  data  set  RN 
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Figure  44.  Correlogram  of  winter  months  only  of 
logged  rainfall  anomalies  from  data  set  RN 
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Figure  45.  Partial  autocorrelations  of  winter 
months  only  logged  rainfall  anomalies  from  data 
set  RN 
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Now,  as  with  the  full  twelve  month  model,  a  look 
at  the  residuals  yields  the  Figures  46  (Residuals) ,  47 
(Correlogram) ,  and  48  (Lag  one  plot) ,  49  (Residuals  vs. 
lag  one) ,  and  Table  19  (Turning  points) .  It  appears  that 
the  residuals  are,  in  fact,  independent.  This  is  similar 
to  the  twelve  month  model. 

The  question  of  the  normality  of  the  residuals  is 
addressed  by  Table  19  and  Figures  50  (Normal  plot)  and  51 
(Histogram) .  The  results  of  these  plots  and  a  basic 
chi-squared  goodness  of  fit  of  22.11  with  17  degrees  of 
freedom  indicate  that  this  winter  month  data  set  is  much 
more  normal  than  was  its  twelve  month  counterpart.  This 
chi-squared  value  is  significant  at  the  .181  level. 
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Figure  46a.  First  order  Markov  residuals  of  logged 
rainfall  anomalies  for  winter  months  only  of  data 
set  RN.  Years  1-12 
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Figure  46b.  First  order  Markov  residuals  of  logged 
rainfall  anomalies  for  winter  months  only  of  data 
set  RN.  Years  12-24 
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Figure  47.  Correlogram  of  first  order  Markov  residuals 
of  logged  rainfall  anomalies  for  winter  months  only  of 
data  set  RN 
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Figure  48.  Lag  one  plot  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  for  winter  months  only  of 
data  set  RN 
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Figure  49.  First  order  Markov  residuals  versus  lag 
one  data  point  from  lagged  rainfall  anomalies  of 
winter  month  only  data  from  data  set  RN 


TABLE  19 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  OF  THE  LOGGED  RAINFALL  ANOMALIES 
OF  DATA  SET  RN 


NUMBER 

OF  TURNING 

POINTS 

■  129 

E[P]  - 

126.667 

V[P] 

=  15. 

PHASE 

LENGTHS 

D 

OBS  . 

EC*] 

1 

32 

78.8 

2 

38 

34.5 

3 

7 

9.9 

4 

1 

2.1 

5 

1 

.4 

6 

0 

0.0 

7 

0 

0.0 
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TABLE  2  0 


GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 
ANOMALIES  OF  DATA  SET  RN 

Moments 


Mean 

Variance 

Skewness 

Kurtosis 


-.001 

.247 

-.161 

-.521 


Percentiles 


Maximum 

-1.148 

Lower  Sixteenth 

-  .845 

Lower  Eight 

-  .616 

Lower  Qurtile 

-  .368 

Median 

.025 

Upper  Quartile 

.  339 

Upper  Eight 

.591 

Upper  Sixteenth 

.  762 

Maximum 

1.  078 
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Figure  50.  Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
winter  months  only  from  data  set  RN 
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FIRST  ORDER  MARKOV  RESIDUAL 5 

Figure  51.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  month  only 
of  data  set  RN 

C.  DATA  SET  FL 

As  in  the  previous  section  on  the  data  sets,  the  analysis 
of  section  III.B  above  carries  forward  fairly  well  to  data 
sets  FL  and  SC.  This  section,  and  the  following,  contain 
only  the  Figures  and  Tables  corresponding  to  those  in  the 
previous  section  on  data  set  RN. 
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1.  Twelve  Month  Series 
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Figure  52.  Partial  correlogram  of  the  logged 
rainfall  anomalies  for  data  set  FL 


TABLE  21 

ESTIMATED  PARTIAL  AUTOCORRELATIONS  FOR 
LOGGED  RAINFALL  ANOMALIES  OF  DATA  SET  FL 


LAG 

VALUE 

LAG 

VALUE 

1 

.185 

14 

-.072 

2 

-.057 

15 

.011 

3 

-.069 

16 

-.019 

4 

.077 

17 

.027 

5 

.016 

18 

-.048 

6 

-.068 

19 

.006 

7 

.010 

20 

-.007 

8 

-.008 

21 

.003 

9 

-.040 

22 

.018 

10 

.022 

23 

-.025 

11 

.031 

24 

.081 

12 

.067 

25 

.012 

13 

.055 

-  • 

Zt  = 

•185Zt-!  +  at 

III. 6 
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Figure  53a.  First  order  Markov  residuals  from 
logged  rainfall  anomalies  of  data  set  FL. 
Months  1-296 
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Figure  53b.  First  order  Markov  residuals  from 
logged  rainfall  anomalies  of  data  set  FL. 
Months  297  -  444, 
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Figure  54.  Autocorrelations  of  residuals  from  first 
order  Markov  process  applied  to  the  logged  rainfall 
anomalies  of  data  set  FL 
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Figure  55.  Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  FL 
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Figure  56.  First  order  Markov  residuals  versus 
lag  one  data  points  from  logged  rainfall  anomalies 
of  data  set  FL 
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TABLE  22 


ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  FROM  DATA  SET  FL 


NUMBER 

OF  TURNING 

POINTS 

=  294 

E[P]  - 

294.667 

V[P] 

=  15. 

PHASE  : 

LENGTHS 

D 

OBS. 

El*] 

1 

188 

183.8 

2 

79 

80.7 

3 

19 

23.2 

4 

6 

5.0 

5 

0 

.9 

6 

1 

.1 

7 

1 

0.0 

8 

0 

0.0 

9 

0 

0.0 

10 

0 

0.0 

TOTALS 

294 

293.7 

TABLE  23 

GENERAL  STATISTICS  OF  FIRST  ORDER  MARKOV 
RESIDUALS  FROM  LOGGED  RAINFALL  ANOMALIES 
OF  DATA  SET  FL 

Moments 


Mean  .000 
Variance  .163 
Skewness  -.044 
Kurtosis  .717 

Percentiles 

Minimum  -1.276 
Lower  Sixteenth  -  . 702 
Lower  Eight  -  .426 
Lower  Quartile  -  .156 
Median  -  .012 
Upper  Quartile  .184 
Upper  Eight  .481 
Upper  Sixteenth  .648 
Maximum  1.124 
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Figure  57.  Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
data  set  FL  _ 


Figure  58.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  FL 


This  data  set  yielded  a  chi-square  value  of  107.66 
for  17  degrees  of  freedom.  The  significance  of  the  value 
is  in  the  zero  plus  range. 

2.  Winter  Series 
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Figure  59a.  Years  1  -  12  of  winter  months  only  of 
logged  rainfall  anomalies  of  data  set  FL 
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Figure  59b.  Years  13  -  37  of  winter  months  only 
logged  rainfall  anomalies  of  data  set  FL 
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Figure  60.  Correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  FL 
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Figure  61.  Partial  correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  FL 

These  displays  indicate  a  model  like 

Z "  -  .199Z^1  +  afc.  III. 7 
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Figure  62a.  Years  1-24,  first  order  Markov 
residuals  of  logged  rainfall  anomalies  for  winter 
months  only,  data  set  FL 
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Figure  62b.  Years  25  -  37,  first  order  Markov 
residuals  of  logged  rainfall  anomalies  for 
winter  months  only,  data  3et  FL 
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Figure  6  3.  Correlogram  of  first  order  Markov 
residuals  of  lagged  rainfall  anomalies  from  winter 
months  only,  data  set  FL 
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Figure  64.  Lag  one  plot  cf  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of  winter 
months  only,  data  set  FL 
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Figure  65.  First  order  Markov  residuals  versus  lag 
one  data  point  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  FL 
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TABLE  24 


ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  OF  THE  LOGGED  RAINFALL  ANOMALIES 
OF  DATA  SET  FL 

NUMBER  OF  TURNING  POINTS  =  209 
E[P]  =  1^6  V[P]  =  15.902 

PHASE  LENGTHS 


D 

OBS  . 

E  [*] 

1 

138 

122.1 

2 

59 

53.5 

3 

8 

15.4 

4 

2 

3.3 

5 

1 

.6 

6 

0 

.1 

7 

0 

0.0 

8 

0 

0.0 

TABLE  25 

GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED  RAINFALL 
ANOMALIES  OF  WINTER  MONTHS  ONLY, 
DATA  SET  FL 


Moments 


Mean 

.000 

Variance 

.234 

Skewness 

-.079 

Kurtosis 

-.376 

Percentiles 

Minimum 

-1.274 

Lower  Sixteenth 

-  .798 

Lower  Eight 

-  .570 

Lower  Quartile 

-  .315 

Median 

.011 

Upper  Quartile 

.  323 

Upper  Eight 

.551 

Upper  Sixteenth 

.748 

Maximum 

1.  136 
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Figure  66.  Standardized  normal  plot  of  first 
order  Markov  residuals  from  logged  rainfall 
anomalies  of  winter  months  only,  data  set  FL 
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Figure  67.  Histogram  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  FL 


The  chi-squared  was  calculated  at  15.35  for  17 
degrees  of  freedom.  This  is  a  significance  level  of  .570, 
thus  indicating  possible  normality. 

D.  DATA  SET  SC 

1.  Twelve  Month  Series 
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Figure  68.  Partial  correlogram  of  the  logged 
anomalies  for  data  set  SC 
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TABLE  26 


ESTIMATED  PARTIAL  AUTOCORRELATIONS 
FOR  LOGGED  RAINFALL  ANOMALIES  OF 
DATA  SET  SC 


LAG 

VALUE 

LAG 

VALUE 

1 

.096 

14 

-.066 

2 

-.075 

15 

.004 

3 

-.053 

16 

-.022 

4 

.046 

17 

.007 

5 

-.004 

18 

-.014 

6 

-.061 

19 

.004 

7 

.042 

20 

-.019 

8 

-.016 

21 

.002 

9 

-.061 

22 

.106 

10 

-.000 

23 

-.046 

11 

-.001 

24 

.059 

12 

.079 

25 

-.015 

13 

.084 

This  information  yields  the  model  as 
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Figure  69a.  First  order  Markov  residuals  from  logged 
rainfall  anomalies  of  data  set  SC.  Months  1  -  148. 
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Figure  69b.  First  order  Markov  residual  from 
logged  rainfall  anomalies  of  data  set  SC. 
Months  149  -  444 
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Figure  69c.  First  order  Markov  residual  from  logged 
rainfall  anomalies  of  data  set  SC.  Months  445  -  596. 
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Figure  70.  Autocorrelations  of  residuals  from  first 
order  Markov  process  applied  to  the  logged  rainfall 
anomalies  of  data  set  SC 
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Figure  71.  Lag  one  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of 
data  set  SC 
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Figure  72.  First  order  Markov  residuals  versus 
lag  one  data  points  from  logged  rainfall  anomalies 
of  data  set  SC 
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TABLE  27 


ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCIES  FROM  THE  FIRST  ORDER  MARKOV 
RESIDUALS  OF  DATA  SET  SC 


NUMBER  OF  TURNING  POINTS 

*  367 

E  [P ]  =*  382.667 

v[p] 

-15.9 

PHASE  LENGTHS 

D 

OBS  . 

E  [*] 

1 

226 

238.8 

2 

95 

104.9 

3 

34 

30.1 

4 

8 

6.6 

5 

2 

1.2 

6 

1 

.2 

7 

1 

0.0 

8 

0 

0.0 

9 

0 

0.0 

10 

0 

0.0 

TOTALS 

36  7 

381.7 

TABLE  28 

GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED 
RAINFALL  ANOMALIES  OF  DATA  SET  SC 

Moments 


Mean 

.000 

Variance 

.209 

Skewness 

.057 

Kurtosis 

.805 

Percentiles 

Minimum 

-1.515 

Lower  Sixteenth 

-  .744 

Lower  Eight 

-  .462 

Lower  Quartile 

-  .205 

Median 

-  .030 

Upper  Quartile 

.207 

Upper  Eight 

.506 

Upper  Sixteenth 

.783 

Maximum 

1.344 
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Figure  73.  Standardized  normal  plot  of  first  order 
Markov  residuals  from  logged  rainfall  anomalies  of 
data  set  SC 
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Figure  74.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  data  set  SC 
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This  data  set  yielded  a  chi-square  value  of  273.95 
for  17  degrees  of  freedom.  This  is  equivalent  to  a  signif¬ 
icance  of  zero  plus. 

2.  Winter  Series 
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Figure  75a.  Years  1  -  12  of  winter  months  only, 
logged  rainfall  anomalies  of  data  set  SC. 
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Figure  75b.  Years  13  -  36  of  winter  months  only, 
logged  rainfall  anomalies  of  data  set  SC 
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Figure  75c.  Years  37  -  48  of  winter  months 
only,  logged  rainfall  anomalies  of  data  set  SC 
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Figure  76.  Correlogram  of  winter  months  only, 
logged  rainfall  anomalies  from  data  set  SC 
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Figure  77.  Partial  correlogram  of  winter  months 
only,  logged  rainfall  anomalies  from  data  set  SC 


This  information  indicates  the  model 
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Figure  78a.  Years  -  24  first  order  Markov  residuals 
of  logged  rainfall  anomalies,  for  winter  months  only 
data  set  SC 
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Figure  79.  Correlogram  of  first  order  Markov  residuals 
of  logged  rainfall  anomalies  from  winter  months  only, 
data  set  SC. 
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Figure  80.  Lag  one  plot  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  months  only, 
data  set  SC 
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Figure  81.  First  order  Markov  residuals  versus  lag 
one  data  point  from  logged  rainfall  anomalies  of 
winter  months  only,  data  set  SC 


TABLE  29 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING  POINTS 
AND  ACTUAL  AND  EXPECTED  PHASE  FREQUENCIES  FROM 
THE  FIRST  ORDER  MARKOV  RESIDUALS  OF  THE 
LOGGED  RAINFALL  ANOMALIES  OF  THE  WINTER  MONTHS 
ONLY,  DATA  SET  SC 


NUMBER  OF  TURNING  POINTS  »  32 
E  [P ]  *  30.667  V[P]  *  15.39 


PHASE  LENGTHS 
D 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

TOTALS 


OBS. 

21 

8 

3 

0 

0 

0 

0 

0 

0 

0 

32 


E[*] 

18.8 

8.1 

2.3 

.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

29.7 
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TABLE  30 


GENERAL  STATISTICS  OF  FIRST  ORDER 
MARKOV  RESIDUALS  FROM  LOGGED 
RAINFALL  ANOMALIES  OF  WINTER  MONTHS 
ONLY,  DATA  SET  SC 


Moments 


Me  an 
Variance 
Skewness 
Kurtosis 


.000 
.  305 
.015 
-.  363 


Percentiles 


Minimum 

-1.513 

Lower  Sixteenth 

-  .872 

Lower  Eight 

-  .663 

Lower  Quartile 

-  .359 

Median 

-  .026 

Upper  Quartile 

.355 

Upper  Eight 

.682 

Upper  Sixteenth 

.882 

Maximum 

1.338 
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Figure  82.  Standardized  normal  plot  of  first  order  Markov 
residuals  from  logged  rainfall  anomalies  of  winter  months 
only,  data  set  SC 
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HISTOGRAM  jF 

FIRST  ORDER  MARKOV  RESIDUALS 

Figure  83.  Histogram  of  first  order  Markov  residuals 
from  logged  rainfall  anomalies  of  winter  months  only, 
data  set  SC 

The  chi-square  was  calculated  at  16.60  for  17 
degrees  of  freedom.  This  is  significant  at  the  .482 
level  thus  indicating  probable  normality. 
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IV.  VALIDATION  OF  FIRST  ORDER  MARKOV  MODELS 

A.  THEORY 

The  general  model  proposed  by  a  first  order  Markov 
process  is,  as  stated  before; 

*t  "  p£t  +  at  1V'1 

where  {a.  }  are  distributed  iid  N(0,  a*) .  To  validate 

L  cL 

this  model,  preferably  independent  data  should  be  subjected 

to  the  model,  and  an  analysis  of  the  residual,  or  forecast  i 

errors,  made. 

As  stated  previously,  years  1975  through  1980  were 
reserved  for  the  purpose  of  validation.  The  method  of 

validation  was  to  use  the  model  to  construct  a  series  of  j 

I 

one  step  ahead  forecasts.  Let  et(l)  be  the  error  in  a  || 

'  1 

forecast  of  time  t+1  from  the  model  at  time  t.  Then  the  j 

minimum  mean  squared  error  forecast  (see  Box  and  Jenkins)  is;  ) 

et(l)  -  Zt  -  PZt_1  IV. 2 

If  the  model  is  correct,  the  sequence  (et(l) }  will  be 
independent  normally  distributed  with  mean  zero  and  variance 
a*  in  the  following  sections  the  models  are  applied  to  the 

A 

reserved  data  sets  (which  may  also  be  found  in  the  appendixes) , 
and  these  forecast  errors  are  calculated.  The  forecast  errors 
are  then  analyzed  to  determine  if 

(1)  The  errors  are  serially  independent 
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(2)  The  errors  are  distributed  as  normal  random 
variables  with  mean  zero,  and  variance  of. 

Since  the  residual  analysis  of  the  twelve  month  model 
already  indicates  a  poor  fit,  the  twelve  month  model  will 
not  be  validated.  Only  the  winter  month  models  will  be 
checked  for  validity. 


B.  DATA  SET  RN 

Figures  84  (Raw  data)  and  85  (Logged  anomalies)  display 
the  reserved  data  set.  The  logged  anomalies  were  formed  by 
removing  the  means  of  the  analyzed  data.  Table  4,  not  the 
means  of  the  logged  reserved  data.  This  was  done  to  remove 
any  bias  from  the  validation. 


CflTfl  IS  “ROM 


Figure  84.  Reserved  rainfall  data  for  data  set  RN 


DATA  13  PROM  -tr-l- 


Figure  85.  Logged  rainfall  anomalies  of  reserved  data  set  RN 

The  forecast  errors.  Figure  86,  their  correlogram. 

Figure  87,  and  independence  tests.  Table  31  indicate  that 
the  errors  are  indeed  independent. 
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Figure  86 .  Forecast  errors  from  first  order  Markov 
model  applied  to  winter  months  of  reserved  data  set  RN 
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Figure  87.  Correlogram  of 
order  Markov  model  applied 
data  set  RN 
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forecast  errors  from  first 
to  winter  months  of  reserved 


TABLE  31 


ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING 
POINTS  AND  ACTUAL  AND  EXPECTED  PHASE 
FREQUENCE IS  FOR  THE  FORECAST  ERRORS  OF 
THE  FIRST  ORDER  MARKOV  MODEL  APPLIED 
TO  THE  WINTER  MONTHS  OF  RESERVED  DATA 
SET  RN 


NUMBER  OF  TURNING  POINTS  =  32 
E[P]  =  30.667  V[P]  =  15.39 

PHASE  LENGTHS 


D 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

TOTALS 


OBS. 

21 

8 

3 

0 

0 

0 

0 

0 

0 

0 

32 


E  [*] 
18.8 
8.1 
2.3 
.5 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
29.7 


The  normality  of  the  forecast  errors  is  addressed  by 
Table  32  (Statistics) ,  Figure  88  (Normal  plot) ,  89  (Histogram) , 
and  a  simple  chi-squared  test.  The  chi-squared  was  calculated 
as  7.82  with  5  degrees  of  freedom  which  is  significant  at  the 
.167  level.  However,  the  normality  of  the  errors  is  somewhat 
questionable  due  to  the  other  displays. 


TABLE  32 


GENERAL  STATISTICS  OF  FORECAST  ERRORS  FROM 
THE  FIRST  ORDER  MARKOV  MODEL  APPLIED  TO  THE 
WINTER  MONTHS  OF  RESERVED  DATA  SET  RN 

Moments 


Mean  .001 
Variance  .259 
Skewness  -.276 
Kurtosis  -.955 

Percentiles 

Minimum  -1.156 
Lower  Sixteenth  -  .791 
Lower  Eight  -  .634 
Lower  Quartile  -  .397 
Median  .069 
Upper  Quartile  .409 
Upper  Eight  .569 
Upper  Sixteenth  .614 
Maximum  .855 


Figure  88.  Standardized  normal  plot  of  forecast  errors 
from  the  first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  RN 
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Figure  89.  Histogram  of  forecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  sec  RN 


C.  DATA  SET  FL 


As  before,  the  similarity  of  results  for  the  different 
data  set  allows  the  analysis  to  be  portrayed  using  the 
displays  only. 
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Figure  90.  Reserved  rainfall  data  for  data  set  FL 
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Figure  91.  Logged  rainfall  anomalies  of 
reserved  data  set  FL. 
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Figure  92.  Forecast  errors  from  the  first  order 
Markov  model  applied  to  the  winter  months  of 
reserved  data  set  FL 
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Figure  93.  Correlogram  of  forecast  errors  from  first 
order  Markov  model  applied  to  the  winter  months  of 
reserved  data  set  FL 


TABLE  33 

ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING  POINTS  AND  ACTUAL 
AND  EXPECTED  PHASE  FREQUENCIES  FROM  THE  FORECAST  ERRORS 
OF  THE  FIRST  ORDER  MARKOV  MODEL  APPLIED  TO  THE  WINTER 
MONTHS  OF  RESERVED  DATA  SET  FL 

NUMBER  OF  TURNING  POINTS  *  34 


E[P]  *  30.667 

PHASE  LENGTHS 
D 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

TOTALS 


V[P]  =  15.396 


OBS. 

24 

8 

2 

0 

0 

0 

0 

0 

0 

0 

34 


E  [*  ] 
18.8 
8.1 
2.3 
.5 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
29.7 
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TABLE  34 


GENERAL  STATISTICS  OF  FORECAST  ERRORS 
FROM  THE  FIRST  ORDER  MARKOV  MODEL  APPLIED 
TO  THE  WINTER  MONTHS  OF  RESERVED  DATA  SET  FL 

Moments 


Mean  -.015 
Variance  .215 
Skewness  -.155 
Kurtosis  -.970 

Percentiles 

Minimum  -1.031 
Lower  Sixteenth  -  .724 
Lower  Eight  -  .549 
Lower  Quartile  -  . 396 
Median  .049 
Upper  Quartile  .369 
Upper  Eight  .479 
Upper  Sixteenth  .531 
Maximum  .904 


Figure  94.  Standardized  normal  plot  of  forecast 
errors  from  the  first  order  Markov  model  applied 
to  the  winter  months  of  reserved  data  set  FL 


HISTOGPAM  :f 
FORECAST  ERRORS 


Figure  95.  Histogram  of  froecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  FL 


The  chi-squared  statistic  was  calculated  as  12.58 
with  7  degrees  of  freedom,  thus  yielding  a  significance 
level  of  0.083.  This  statistic  and  the  displays  imply 
that  the  data  are  only  marginally  normal  if  at  all. 
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DATA  SET  SC 


Figure  96.  Reserved  rainfall  for  data  set  Sc 


set  SC 
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Figure  98.  Forecast  errors  from  first  order 
Markov  model  applied  to  the  winter  months  of 
reserved  data  set  SC 
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Figure  99.  Correlogram  of  forecast  errors 
from  first  order  Markov  model  applied  to  the 
winter  months  of  reserved  data  set  SC 


TABLE  35 


ACTUAL  AND  EXPECTED  NUMBER  OF  TURNING  POINTS  AND 
ACTUAL  AND  EXPECTED  PHASE  FREQUENCES  FROM  THE 
FORECAST  ERRORS  OF  THE  FIRST  ORDER  MARKOV  MODEL 
APPLIED  TO  THE  WINTER  MONTHS  OF  RESERVED  DATA  SET  SC 


NUMBER  OF  TURNING  POINTS 
E  [P ]  =  30.667 
PHASE  LENGTHS 


34 


D 

1 

2 

3 

4 

5 

6 
7 

TOTALS 


OBS. 

24 

8 

2 

0 

0 

0 

0 

34 


V[P] 


EE*] 

18.8 

8.1 

2.3 

.5 

0.0 

0.0 

0.0 

29.7 


15.396 


TABLE  36 

GENERAL  STATISTICS  OF  FORECAST  ERRORS 
FROM  THE  FIRST  ORDER  MARKOV  MODEL  APPLIED 
TO  THE  WINTER  MONTHS  OF  RESERVED  DATA  SET  SC 


Momen  ts 


Mean  -.003 
Variance  .296 
Skewness  -.298 
Kurtosis  -.413 

Percentiles 

Minimum  -1.386 
Lower  Sixteenth  -  .846 
Lower  Eight  -  .559 
Lower  Quartile  -  .434 
Median  .031 
Upper  Quartile  .414 
Upper  Eight  .588 
Upper  Sixteenth  .735 
Maximum  1.005 
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Figure  100.  Standardized  normal  plot  of  forecast 
errors  from  the  first  order  Markov  model  applied  to 
the  winter  months  of  reserved  data  set  SC 


Figure  101.  Histogram  of  forecast  errors  from  the 
first  order  Markov  model  applied  to  the  winter 
months  of  reserved  data  set  SC 
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The  chi-squared  statistic  was  calculated  as  9.71  with. 
7  degrees  of  freedom,  thus  yielding  a  significance  level 
of  .205. 


E.  CONCLUSIONS 

The  application  of  a  Markovian  model  was  indicated  by 
the  apparent  dependence  of  adjacent  months  and  the  apparent 
lack  of  dependence  at  any  other  lag.  The  preceeding  sub¬ 
sections,  however,  indicate  that  the  first  order  Markovian 
model  is  weak  at  best. 

The  structure  of  the  data,  visually,  still  points 
toward  some  sort  of  underlying  order.  The  following 
sections  attempt  to  discover  this  order. 
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V.  2x2  TABLES 


A.  THEORY 

As  seen  in  sections  III  and  IV,  the  classical  ARMA  time 
series  approach  does  not  seem  to  adequately  describe  the 
data.  Another  technique  used  to  explore  possible  relation¬ 
ships  is  the  2x2  contingency  tables. 

The  idea  to  be  explored  is  whether  or  not  some  subset  of 
the  data,  to  be  called  the  control,  may  be  used  to  predict 
in  some  way  the  behavior  of  another  subset  of  the  data,  to 
be  called  the  complement.  Here,  the  data  are  reduced  from 
monthly  observations  to  yearly  observations  as  described 
below. 

I 

Let  X  be  the  subset  of  a  year,  to  be  called  the  control, 
and  let  Y  be  the  subset  to  be  called  the  complement.  It  is 

—  i  i 

necessary  that  XOY  =  0;  that  is,  the  intersection  of  these 

) 

two  sets  is  empty.  The  data  are  then  compared  for  some 

quality  in  X  and  for  some  quality  in  Y.  The  question  is  I 

then:  does  the  presence  (or  absence)  of  the  quality  in  X  j 

affect  the  presence  (or  absence)  of  the  quality  in  Y?  An  j 

example  of  a  typical  table,  is  shown  below  in  Figure  102.  j 
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•  • 

Figure  102. 


Typical  2x2  contingency  table 


The  table  elements,  n^,  represents  the  number  of  years 
which  display  quality  i  in  the  control  and  quality  j  in 
the  complement.  The  marginal  entries  n^.  and  n.^  represents 
the  numbers  of  years  for  which  the  control  has  quality  i  and 
the  number  of  years  the  complement  has  quality  j  respectively. 
The  overall  number  of  years,  n. . ,  is  in  the  lower  right  of 
the  table. 

Brownlee  [Ref.  6]  contains  a  very  good  discussion  of  the 
theory  and  use  of  2x2  contingency  tables.  Using  the  notation 
of  Brownlee,  let  9^  be  the  probability  that  any  given  year 
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will  have  a  control  quality  i  and  a  complement  quality  j. 
Then  estimates  of  the  Q^'s  are 


A 


“ij/”- 
n ,  /n . . 

n.  j/n 


V.l 


If  the  control  and  complement  are  independent. 


V.  2 


These  simple  assumptions  allow  for  a  thorough  investigation 
of  the  possible  interrelationships  within  the  data  sets. 

Another  way  to  view  the  assumption  of  independence  is 
through  the  use  of  proportions.  Thus,  if  the  basic  division 
is  made  via  the  quality  of  the  control,  the  the  proportions 
*  nn//ni*  (respectively  P2  *  n21^n2**  V.l 

represent,  in  words,  the  proportion  of  the  years  that  have 
quality  1  in  the  control  and  have  quality  1  (respectively 
2)  in  the  complement. 

The  question  of  independence  may  be  approached  in  several 
ways  as  described  below. 

1.  Fishers  Exact  Test 

A  test  for  the  significance  of  any  dependence  was 
proposed  by  Fisher  in  the  case  in  which  the  marginal  totals, 
ni*»n,j,  n. .  are  known  a  priori  (cf.  Brownlee  [Ref.  6]). 
To  draw  from  Brownlee,  knowledge  of  the  marginal.'  and  n^ 
gives  knowledge  of  all  the  other  elements  of  the  table.  The 
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probability  of  the  event  of  having  exactly  n^  years  that 
display  quality  1  in  both  the  complement  and  the  control  is; 


p(NU=nll> 


n, 


inllinl2 ln21!n22! 


V.  4 


A  test  may  then  be  applied,  using  V.4,  to  determine  the 
significance  of  any  dependence.  This  test  is  usually  applied 
by  simply  summing  these  probabilities  in  the  tail  of  the 
distribution  (V.4)  in  the  same  direction  as  the  noted  extreme. 

The  usual  procedure  to  provide  a  two-sided  test  of 
significance  is  to  double  a  one  sided  figure.  This  procedure 
is  acceptable  due  to  the  symmetric  appearance  of  the  distri¬ 
bution. 

Under  the  assumptions  of  independence 


and  the  random  variable 


N  -E[N 

U  * 

'/V[N11] 


U 

11 


] 


defined  as 


V.  7 


is  asymptotically  distributed  as  a  normal  random  variable 
with  mean  zero  and  variance  one.  This  asymptotic  result 
combined  with  a  continuity  correction  yields  a  test  statistic 
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of 


The  statistic  U'  may  then  be  used  as  a  test,  using 
standard  normal  tables,  of  the  significance  of  any  variation 

from  the  assumption  of  independence.  It  should  be  noted  at 

2 

this  point  that  if  the  random  variable  U'  is  squared,  tJ' 
will  be  distributed  as  a  chi-square  with  one  degree  of  freedom. 
The  squaring  of  U*  with  simplifying  algebra  yields  the  Yates 
correction  to  a  standard  chi-squared  goodness  of  fit  statistic 


(u‘)2 


{  I ^22  ™^12  ^21  I •  •  /2  }  n.  . 

(nH+nl2)  (nn+n2i)  (ni2+n22)  {n21+n22‘ 


see  Dixon  and  Massey  [Ref.  5],  This  allows  the  use  of  the 
chi-square  tables  as  an  equivalent  test  to  that  of  V.8. 

2 .  Odds 

Subsection  V.A.l  above  deals  with  the  significance 
of  any  observed  interdepedence  between  the  control  and  the 
complement.  The  question  of  the  degree  of  dependence  should 
also  be  addressed.  The  measure  to  be  used  is  the  odds  ratio. 


Using  the  notation  of  Fleiss  [Ref.  3],  a  measure  of  seeing 
quality  1  in  the  complement  Y  may  be 


-  P(Yall  X«l) 
1  "  P (Y=2  X-l) 


V.10 


this  is  then  the  odds  that  quality  1  will  occur  in  the  comple¬ 
ment  given  that  quality  1  is  present  in  the  control.  In  a 
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B.  ANALYSIS 


The  theory  of  subsection  A  above  is  applied  to  the  three 
data  sets  as  discussed  below.  The  control  is  typically  taken 
as  a  monthly  anomaly,  say  October.  Here,  the  quality  is  taken 
as  either  a  positive  or  a  negative  anomaly.  Thus  X=1  occurs 
when  the  month  of  October  falls  below  its  mean  and  X=2  occurs 
when  it  falls  above  its  mean.  The  complement  consists  of  the 
sum  of  the  rainfall  for  the  succeeding  eleven  months,  or  in 
symbols : 


12 

*t-  l 


m=2 


t  ,ra 


1 

N 


N 

l  ( 


t=l 


12 

l  _) 

m=2  t,m 


V.  15 


Where  it  is  understood  that  X=1  when  Xfc<  0,  X=2  when 

Xfc  >  0  and  similarly  for  Yfc. 

Various  control  subsets  are  used;  October  through 
September  were  investigated  by  themselves  as  were  all 
adjacent  pairs,  triples,  and  four-tuples  of  months.  For 
an  example,  consider  the  spring  (April,  May,  and  June)  and 
its  complement  (July  through  March) .  In  this  case 


N-l 


Xt  =  £  Rt  m  "  i  l  (  l  R  ) 

m»7  N  t-1  'm*7  Rt,m' 


V.  16 


12 


N-l  12 


’  m-10Rt'm  +  Jl  ^1-  *  ^  Jl'JlO1'”  + 

Equations  V.  15  and  V.16  imply  that  the  data  are  always 
analyzed  as  deviations  from  the  arithmetic  mean.  However, 
the  data  are  also  analyzed  as  deviations  from  the  median 
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and  the  lower  quartile.  In  the  tables  to  follow,  ’A'  refers 
to  both  control  and  complement  having  the  arithmetic  mean 
removed,  'M'  refers  to  both,  control  and  complement  having 
their  respective  medians  removed,  and  'QL'  refers  to  the 
control  having  the  lower  quartile  removed  while  the  median 
was  removed  from  the  complement. 

The  first  four  Tables  (37  through  40) ,  give  the  signifi¬ 
cance  levels  of  observed  departures  from  independence  of  the 
control  and  complement.  Only  those  values  having  a  Yates 
corrected  chi-square  of  greater  than  1.00  are  listed.  The 
entries  represent  the  two-tailed  probability  of  a  random 
deviation  in  excess  of  that  observed.  Although  the  cut  off 
criterion  was  the  Yates  chi-square,  the  agreement  between 
its  probability  and  that  obtained  from  the  Fisher  exact  and 
normal  tests  did  not  differ  in  the  first  two  decimal  places. 


TABLE  37 

SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM  INDEPENDENCE  OF 
SINGLE  MONTH  CONTROL  VERSUS  SUCCEEDING  ELEVEN 
MONTH  COMPLEMENTS 


Data  set 

RN 

FL 

SC 

Differentiator  A 

M 

QL 

A 

M 

QL 

A 

M 

Control 

October  . 18 

.22 

.24 

.31 

in 

H 

• 

November 

.10 

.12 

December 

January  . 14 

.30 

.19 

.002 

.04 

.12 

.21 

.19 

February 

.11 

March 

.13 

April 

.18 

May 

.14 

June  .24 

July 

.31 

.21 

August 

September 

.14 

.12 
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TABLE  38 


SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM  INDEPENDENCE 
OF  PAIRS  OF  MONTHS  VERSUS  SUCCEEDING 
TEN  MONTH  COMPLEMENTS 


Data  Set 

RN 

FL 

SC 

Differentiators 

A 

M 

QL 

A 

M 

QL 

A 

M 

QL 

Control 

Oct+Nov 

.18 

.09 

.24 

Nov+Dec 

.10 

.12 

.15 

.  31 

.08 

Jan+Feb 

Feb+Mar 

.30 

.  30 

.01 

.02 

.12 

.  11 

.19 

.11 

Mar+Apr 

Apr+May 

May+Jun 

Jun+Jul 

.14 

.12 

.24 

Jul+Aug 

Aug+Sep 

.19 

00 

H 

• 

TABLE  39 

SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM  INDEPENDENCE 
OF  TRIPLES  OF  MONTHS  VERSUS  SUCCEEDING 
NINE-MONTH  COMPLEMENTS 

Data  Set  RN  FL  SC 

Differentiators  A  M  QL  A  M  QL  A  M  QL 

Control 
Oct+Nov+Dec 
Nov+Dec+Jan 
Dec+Jan+Feb 
Jan+Feb+Mar 
Feb+Mar+Apr 
Mar+Apr+May 
Apr+May+Jun 
May+Jun+Jul 
Jun+Jul+Aug 
Jul+Aug+Sep 


19 


.31 


.31  .28 


.22 


.30 

.14  .30 


.18 
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TABLE  40 


SIGNIFICANCE  OF  OBSERVED  DEPARTURES  FROM 
INDEPENDENCE  OF  FOUR-TUPLES  OF  MONTHS  VERSUS 
SUCCEEDING  EIGHT-MONTH  COMPLEMENTS 


Data  set 

RN 

FL 

SC 

Differentiator 

A  M 

QL 

A  M  QL  A 

M 

QL 

Control 

Oct+Nov+Dec+Jan 

Nov+Dec+Jan+Feb 

.22 

.26 

Dec+Jan+Feb+Mar 

Jan+Feb+Mar+Apr 

.22 

.31 

Feb+Mar+Apr+May 

.22 

.27 

Mar+Apr+May+Jun 

Apr+May+Jun+Jul 

May+Jun+Jul+Aug 

Jun+Jul+Aug+Sep 

.14 

.30 

.20 

.  31 

.28 

Several  choices  for  predictors  are  suggested  in  the 
previous  tables.  However,  the  apparent  strongest  candidate 
for  a  predictor  is  January.  The  control  of  January  by 
itself  and  January  paired  with  December,  are  the  most 
consistently  significant  entries.  Tables  41  below  gives 
the  odds  ratio,  V.13,  for  January,  and  January  and  December, 
as  controls. 
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TABLE  41 


ODDS  RATIO  OF  JANUARY  VERSUS  FEBRUARY 
THROUGH  DECEMBER  AND  JANUARY  PLUS 
DECEMBER  VERSUS  FEBRUARY  THROUGH  NOVEMBER 

Differentiator  A  M  QL 

Data  set  RN 

January  4.59  3.15  5.13 

Jan+Dec  3.15  3.15  2.01 

Data  set  FL 

January  10.71  4.72  4.30 

Jan+Dec  7.42  6.02  4.30 

Data  set  SC 

January  2.45  2.49  2.23 

Jan+Dec  3.00  2.49  1.44 

At  this  point  in  the  analysis  it  was  decided  to  explore 
more  fully  the  power  of  January  as  a  predictor.  It  should 
be  stated  that  other  possibilities  for  predictors  are 
suggested  by  the  tables,  but  time  did  not  allow  am  exhaustive 
study  of  all  of  these. 

C.  OTHER  RESULTS 

The  results  of  section  V.B  suggest  that  a  more  detailed 
analysis  of  January  as  a  predictor  is  in  order.  The  first 
method  tried  for  this  was  ordinary  least  squares  regression 
of  the  rainfall  total  in  January  versus  the  total  for  February 
through  December.  This  is  the  model  below. 

Let  Xt  =  R  4 


then  assume  that 


Yfc  =  a  +  8Xfc  +  £ t  V. 18 

as  the  standard,  linear  model  where  {£fc}  are  assumed  to  be 
independent  and  identically  distributed  with  mean  zero  and 
variance  ae2.  If  the  predictability  of  January  is  strong, 
this  model,  V.18,  may  result  in  a  good  fit  of  the  data. 

Table  42  below  is  the  resulting  ANOVA  for  this  regression. 

As  may  easily  be  seen,  the  model  does  not  appear  to  have  any 
significance. 


TABLE  42 

ANOVA  FOR  REGRESSION  OF  SIMPLE 
LINEAR  MODEL  FOR  ALL  DATA  SETS 

R-squared  =  .1033 

Standard  error  of  estimate  »  4.1574  RN 


SOURCE 

DF 

SS 

MS 

F 

Total 

22 

404.798 

Regression 

1 

41. 821 

41.821 

2.42 

Jan -Control 

1 

41. 821 

41.821 

2.42 

Residual 

21 

362.977 

17.285 

Variable 

Coefficient 

Standard  error 

T 

Alpha 

.2178 

.878 

.25 

Beta 

.5993 

.  385 

1.56 

90%  Confidence 

Limits 

Lower  limit 

Upper 

■  limit 

Alpha 

-1.1817 

1. 

6174 

Beta 

-  .0148 

1. 

2133 
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R-squared  *  .2066 

Standard  error  of  estimate  »  4.28Q4 


FL 


AOV 

SOURCE 

DF 

SS 

MS 

F 

Total 

35 

785.196 

Regression 

1 

162.252 

162.252 

8.86 

Jan -Control 

1 

162.252 

162.252 

8.86 

Residual 

34 

622.945 

18.322 

Variable 

Coefficient 

Standard  error  T 

Alpha 

-.2670 

.715 

-.37 

Beta 

1.0095 

.  339 

2.98 

90%  Confidence 

Limits 

Lower  limit 

Upper  limit 

Alpha 

-  .2670 

.8543 

Beta 

1.0095 

1.5417 

R-squared  = 

.0426 

Standard  error  of  estimate  = 

6. 1508 

SC 

AOV 

SOURCE 

DF 

SS 

MS 

F 

Total 

46 

1778.131 

Regression 

1 

75,687 

75.687 

2.00 

Jan-Control 

1 

75.687 

75.687 

2.00 

Residual 

45 

1702.444 

37.832 

Variable 

Coefficient 

Standard  error  T 

Alpha 

-.2771 

.898 

-.31 

Beta 

.4355 

.  308 

1.41 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

-1.6766 

1.1224 

Beta 

-  .0446 

.9157 

The  same 

model  with  means  removed,  V 

.19,  below  was  tried 

and,  although  slightly  better,  is  still 

not  strong. 

- 

ct  +  6(Xt-X.)  + 

£t 

V.  19 
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TABLE  43 


ANOVA  FOR  REGRESSION  OF  SIMPLE  LINEAR 
MODEL  WITH  MEANS  REMOVED  FOR  ALL  DATA  SETS 


R-squared  = 

.1033 

Standard  error  of  estimate 

=  4.1575 

AOV 

SOURCE 

DF 

SS 

MS 

Total 

22 

404.798 

Regression 

1 

41.821 

41.821 

Jam -Control 

1 

41.821 

41.821 

Residual 

21 

362.977 

17.285 

Variable 

Coefficient 

Standard  error 

Alpha 

11.7215 

1.851 

Beta 

.5993 

.  385 

90%  Confidence  Limits 

Lower  limit 

Upper 

limit 

Alpha 

8.7236 

13. 

2483 

-.0148 

1. 

2133 

R-squared  = 

.1747 

Standard  error  of  estimate 

*  4.4276 

AOV 

SOURCE 

DF 

SS 

MS 

Total 

35 

807.618 

Regression 

1 

141.108 

141.108 

Jan -Control 

1 

141.108 

141.108 

Residual 

34 

666.510 

19.603 

Variable 

Coefficient 

Standard  error 

Alpha 

10.9860 

1. 

442 

Beta 

.9414 

• 

351 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

8. 7236 

13.2483 

Beta 

.  3909 

1.4920 

RN 


F 


2.42 

2.42 


T 

6.33 

1.56 


FL 


F 


7.20 

7.20 


T 

7.62 

2.68 
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R-squared  =  .5617  SC 

Standard  error  of  estimate  =  .7589 


AOV 

SOURCE 

DF 

ss 

MS 

F 

Total 

7 

7.884 

Regression 

1 

4.428 

4.428 

7.69 

Jan -Control 

1 

4.428 

4.428 

7.69 

Residual 

6 

3.456 

.576 

Variable 

Coefficient 

Standard  error 

T 

Alpha 

-  .1946 

.287 

-.68 

Beta 

1.7710 

.637 

2.77 

90%  Confidence  : 

Limits 

Lower  limit 

Upper  limit 

Alpha 

-.7040 

.  3148 

Beta 

.6362 

2.9058 

The  amount  of  rainfall  in  January  does  not  appear  to  be 
a  strong  predictor  for  the  amount  of  rainfall  in  February 
through  December.  This  seems  to  indicate  that  the  relation¬ 
ships  between  January  rainfall  and  rainfall  during  the  next 
eleven  months  is  not  as  strong  as  expected.  However,  a 
further  technique  is  available,  that  of  log-odds  and  logistic 
regression,  which  are  the  subjects  of  the  next  section. 
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VI.  LOGISTIC  ANALYSIS 


A.  THEORY 


The  logistic  analysis  to  be  described  in  the  section  was 
developed  from  Gaver  [Ref.  2]  and  Fleiss  [Ref.  3],  This 
analysis  derives  from  the  model  of  1.5  as  stated  in  the 
introduction . 

The  basic  approach  is  to  view  the  complement  as  having  a 
binary  representation,  with  success  being  defined  as  a  comple¬ 
ment  above  its  mean  (see  equations  1.3  and  1.4)  and  failure 
as  the  complement  below  its  mean.  The  problem  then  is  to  find 
the  conditional  probability  of  a  success  (the  complement  being 
above  its  mean  for  a  year)  given  that  the  control  (January 
rainfall)  takes  on  a  particular  value. 

The  control  is  now  taken  to  be  the  logged  rainfall  anomaly 
of  January  and  is  found  in  equation  1.2  repeated  below. 


X 


t 


N 

l  4) 

t=l 


VI.  1 


If  the  probability  of  success,  given  Xfc  is  written  as 

P (Success j Xfc)  =  0  VI.  2 

a  superficially  attractive  model  for  0fc  is 

9t  *  a  +3  Xt  +  £ t  VI. 3 

This  model  has  two  difficulties,  the  worst  of  which  is 
that  probabilities  of  greater  than  one  or  less  that  zero  are 
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allowed.  Secondly,  the  8^  are  not  available  in  proportion 
form  to  fit  the  model. 

Initially,  the  problem  of  estimating  is  approached 

by  grouping  the  data.  Section  II  indicated  that  each  year 
seemed  to  be  independent  of  the  next  and  that  there  was  no 
trend.  This  allows  for  the  ordering  of  the  Xt's  into  their 
order  statistics  X^'s.  Once  the  ordering  has  been  done, 
non  overlapping  groups  of  arbitrary  size  may  be  formed  as 
shown.  Let  X^,X2,X2» . . • ,  X^  a  series  °f  12  years  from  an 
arbitrary  data  set,  with  associated  order  astatistics 
X ( i ) »X(2) ,x(3) • •• • ,x(i2) *  T^e  ^  groups  of  size  three  are 
desired,  the  data  are  partitioned  below. 


X(l) •  X (2) ,X( 3) lX(4) ,X(5) ,X(6) 


...  X 


(10) ,x(ll) ,X(12) 


Given  these  groups,  let  X..  be  a  measure  of  location  for 
the  jth  group.  This  analysis  used  the  median,  therefore 
Xj  =  X(3j+1)’  Also  associated  with  each  group  is  Rj  (not 
rainfall),  the  number  of  success  in  group  j,  and  n ^ ,  the 
number  of  elements  in  group  j.  From  this  set  up,  the  required 
probabilities  may  be  estimated  as; 

VI.  4 

A  solution  to  the  first  problem,  that  of  the  model  yield¬ 
ing  probabilities  outside  of  (0,1),  is  to  use  the  log  odds, 
instead  of  9j  where; 


9 .  =  R./n . 
3  3  3 


Log  odds 


l-Tjy 


VI. 5 
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1 


which,  is  equivalent  to  the  logarithm  of  the  odds  ratio  as 
given  in  V.13.  Gaver  [Ref.  2]  suggests  that  a  correction  of 
.5  be  applied  to  guard  against  the  problem  of  0  and  1  within 
the  logarithm  and  to  reduce  the  bias.  The  statistic  then 
becomes; 

*  in 

The  temptation  is  to  go  directly  to  the  model 

4>  j  *  ct  +  3  Xj ;  VI.  7 

yet  V[9j]  =  0j(l~6j)/nj  which  is  not  constant.  This 
suggests  the  need  for  a  weighting  scheme. 

The  weighting  scheme  used  was  that  of  iteratively 
reweighted  least  squares,  (IRWLS) ,  using  the  bi-weights. 

This  method  is  dicussed  in  detail  in  Mosteller  and  Tukey 
[Ref.  13]. 

Although  grouping  of  the  data  and  the  model  of  VI. 7 
provide  adequate  representation  of  the  underlying  structure, 
the  logistic  model  itself,  1.5,  when  viewed  through  the  eyes 
of  maximum  likelihood  theory  can  yield  more  insight. 

The  model  is  assumed  to  be 

ec*8  Xt 

3C  =  —  ea-BXt  VI. 3 

where  the  Xt  are  independent. 


VI.  6 
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The  likelihood  function  is  then 


N 

L  ( X ,  Y ;  a ,  3 )  =  J~[ 
t=l 


rr  e«V6xtYt 

t*l  n  (l+ea+SXt) 


and  the  log-likelihood  is 
M  N 

L(X, Y;a,0)  -  a  %  Y.  +  s£  X.Y. 

t-1  t-1  C 

-  !  In  (1  +  e“t6xt) 
t-1 

The  gradient,  and  Hessian  of  VI. 9  are 
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A  simple  way  to  solve  for  a  and  3  is  to  use  Newtons 
Method  at-. 


CO 

e/k+i 


H^1  V  tL 


since  all  necessary  elements  may  be  calculated  in  one  pass  of 
the  computor  algorithm. 

One  beneficial  byproduct  of  the  maximum  likelihood 
approach  is  the  asymptotic  information  matrix,  H^.  Gaver 
[Ref.  2]  states  that  the  diagonal  elements  of  this  matrix 
provide  good  estimates  of  V[aj  and  V[8]  under  assumptions 
of  normality. 


B.  ANALYSIS 

1.  Grouped  Data 

The  first  approach  taken  was  that  of  grouping  the 
data  as  described  above.  Groups  of  3,  4,  and  5  were  used, 
as  were  two  separate  methods  of  regression,  ordinary  least 
squares  (OLS )  ,  and  iteratively  reweighted  least  squares 
( IRWLS) .  Tables  44  (RN) ,  45  (FL) ,  and  46  (SC)  present  the 
data  and  Tables  47  (OLS) ,  and  48  (IRWLS)  present  the  results 
of  the  regressions. 
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TABLE  44 


DATA  SET  RN,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  FORMS 


UNGROUPED 

GROUPED 

t 

X 

Y 

j 

X 

Y 

1 

-1.22 

0 

group® 3 

2 

-1.09 

0 

1 

-1.09 

0 

3 

-  .66 

0 

2 

-  .46 

1 

4 

-  .58 

0 

3 

-  .19 

1 

5 

-  .46 

0 

4 

.01 

2 

6 

-  .36 

1 

5 

.16 

2 

7 

-  .31 

0 

6 

.36 

2 

8 

-  .18 

1 

7 

.50 

2 

9 

-  .17 

0 

8 

.98 

2 

10 

-  .05 

1 

group® 4 

11 

.01 

1 

1 

-.  38 

0 

12 

.01 

0 

2 

-.34 

2 

13 

.15 

1 

3 

-.03 

2 

14 

.16 

0 

4 

.20 

2 

15 

.24 

1 

5 

.47 

4 

16 

.28 

0 

6 

.94 

2 

17 

.36 

1 

group® 5 

18 

.46 

1 

1 

-.66 

0 

19 

.48 

1 

2 

-.19 

3 

20 

.50 

1 

3 

.15 

3 

21 

.51 

0 

4 

.46 

4 

22 

.94 

1 

5 

.94 

2 

23 

1.01 
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TABLE  45 


DATA  SET  FL,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  DATA 


UNGROUPED 

GROUPED 

t 

X 

Y 

j 

X 

Y 

1 

-3.41 

1 

group= 3 

2 

-1.82 

0 

1 

-1.  82 

1 

3 

-1.14 

0 

2 

-1.00 

0 

4 

-1.02 

0 

3 

-  .46 

0 

5 

-1.00 

0 

4 

-  .20 

2 

6 

-  .92 

0 

5 

.11 

1 

7 

-  .47 

0 

6 

.14 

1 

8 

-  .46 

0 

7 

.20 

0 

9 

-  .39 

0 

8 

.34 

2 

10 

-  .28 

1 

9 

.50 

2 

11 

-  .20 

0 

10 

.61 

2 

12 

-  .14 

1 

11 

.75 

3 

13 

-  .04 

1 

12 

1.14 

3 

14 

.11 

0 

group=4 

15 

.11 

0 

1 

-1.48 

1 

16 

.13 

0 

2 

-  .70 

0 

17 

.14 

0 

3 

-  .24 

2 

18 

.16 

1 

4 

.11 

1 

19 

.16 

0 

5 

.16 

1 

20 

.20 

0 

6 

.32 

2 

21 

.25 

0 

7 

.50 

2 

22 

.30 

1 

8 

.71 

4 

23 

.34 

1 

9 

1.11 

4 

24 

.38 

0 

groups 5 

25 

.46 

1 

1 

-1.14 

1 

26 

.50 

0 

2 

-  .46 

1 

27 

.51 

1 

3 

-  .04 

2 

28 

.57 

0 

4 

.16 

1 

29 

.61 

1 

5 

.  34 

3 

30 

.71 

1 

6 

.57 

3 

31 

.71 

1 

7 

.95 

6 

32 

.75 

1 

33 

.82 

1 

34 

1.08 

1 

35 

1.14 

1 

36 

1.31 

1 
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TABLE  46 


DATA  SET  SC,  LOGGED  JANUARY  ANOMALIES  AND 
SUCCESSES  FOR  GROUPED  AND  UNGROUPED  FORMS 


UN GROUPED 


X  Y 

-4.24  1 

-1.41  0 

-1.38  0 

-1.08  0 

-1.07  0 

-1.06  0 

-1.02  0 

-  .96  1 

-  .65  0 

-  .60  0 

-  .55  1 

-  .36  0 

-  .27  1 

-  .26  0 

-  .21  1 

-  .17  1 

-  .10  0 

-  .03  0 

-  .01  1 

.02  0 

.12  1 

.14  0 

.20  0 

.22  0 

.23  0 

.28  1 

.28  1 

.40  1 

.40  1 

.41  0 

.43  0 

.46  1 

.49  0 

.56  1 

.56  1 

.60  0 

.64  0 

.69  0 

.70  0 

.70  1 

.  78  1 

.81  0 

.  88  0 

.97  0 

.97  1 

1.03  1 

1.47  1 


GROUPED 

j 

X 

Y 

group* 3 

1 

-1.41 

1 

2 

-1.07 

0 

3 

-  .96 

1 

4 

-  .55 

1 

5 

-  .26 

2 

6 

-  .10 

1 

7 

.02 

2 

8 

.20 

0 

9 

.28 

2 

10 

.40 

2 

11 

.46 

1 

12 

.56 

2 

13 

.69 

0 

14 

.78 

2 

15 

.97 

1 

16 

1.25 

2 

group=4 

1 

-1.39 

1 

2 

-1.04 

1 

3 

-  .58 

1 

4 

.23 

3 

5 

-  .02 

1 

6 

.17 

1 

7 

.28 

3 

8 

.42 

2 

9 

.56 

2 

10 

.69 

1 

11 

.78 

1 

12 

1.03 

3 

group* 5 

1 

-1.38 

1 

2 

-  .96 

1 

3 

-  .27 

3 

4 

-  .04 

2 

5 

.20 

1 

6 

.40 

4 

7 

.49 

3 

8 

.69 

1 

9 

.88 

2 

10 

1.25 

2 
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TABLE  47a 


ORDINARY  LEAST  SQUARES  REGRESSION  WITH 
THE  MODEL  OF  EQUATION  VI. 7  FOR  DATA  SET  RN 

R-squared  =  .0423 

Standard  error  of  estimate  =  6.1508  RN 

GROUP* 3 


AOV 


SOURCE 

DF 

SS  MS 

F 

Total 

46 

1778.131 

Regression 

1 

75.687  75.687 

2.00 

Jan -Control 

1 

75.687  75.687 

2.00 

Residual 

45 

1702.444  37.832 

Variable 

Coefficient 

Standard  error 

T 

Alpha 

14.5833 

1.677 

8.69 

Beta 

.4355 

.  308 

1.41 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

11.9680 

17.1988 

Beta 

-  .0446 

.9157 

R-squared  * 

.4068 

RN 

Standard  error  of  estimate 

*  1.2099 

GROUP* 4 

AOV 

SOURCE 

DF 

SS  MS 

F 

Total 

5 

9.871 

Regression 

1 

4.016  4.016 

2.74 

Jan-Control 

1 

4.016  4.016 

2.74 

Residual 

4 

5.855  1.464 

Variable 

Coefficient 

Standard  error 

T 

Alpha 

-.1696 

.517 

-.33 

Beta 

1.7687 

1.068 

1.66 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

-1.1675 

.8284 

Beta 

1.7688 

3.8288 

R-squared  *  .5973 

Standard  error  of  estimate  *  .9995 


RN 

GROUP* 5 


AOV 


SOURCE 

DF 

ss 

MS 

F 

Total 

4 

7.442 

Regression 

1 

4.446 

4.446 

4.45 

Jan -Control 

1 

4.446 

4.446 

4.45 

Residual 

3 

2.997 

.999 

Variable 

Coefficient 

Standard  error 

T 

Alpha 

-.2647 

.461 

-.57 

Beta 

1.7193 

.815 

2.11 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

-1.2358 

.7064 

Beta 

.0040 

3.4345 

TABLE 

47b 

ORDINARY  LEAST 

SQUARES  REGRESSION 

WITH 

THE 

MODEL  OF  EQUATION 

VI.  7 

FOR  DATA 

SET  FL 

R-squared  * 

.2862 

FL 

Standard  error  of  estimate  *  .8949 

GROUP* 

AOV 

SOURCE 

DF 

SS 

MS 

F 

Total 

11 

20.591 

Regression 

1 

7.998 

7.998 

6.35 

Jan-Control 

1 

7.998 

7.998 

6.35 

Residual 

10 

12.593 

1.259 

Variable 

Coefficient 

Standard  error  T 

Alpha 

-.1459 

.324 

-.33 

Beta 

1.0521 

.417 

2.52 

90%  Confidence  Limits 
Lower  limit  Upper  limit 

Alpha  -.6871  .4944 

Beta  .3550  1.7492 
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TABLE  47c 


ORDINARY  LEAST  SQUARES  REGRESSION  WITH 
THE  MODEL  OF  EQUATION  VI . 7  FOR  DATA  SET  SC 


R-squared  =  .1404  SC 

Standard  error  of  estimate  -  .9901  GROUP=3 

AOV 


SOURCE 

DF 

SS 

MS 

F 

Total 

15 

15.965 

Regression 

1 

2.241 

2.241 

2.29 

Jam -Control 

1 

2.241 

2.241 

2.29 

Residual 

14 

13.724 

.980 

Variable 

Coefficient 

Standard  Error 

T 

Alpha 

-.  3039 

.249 

-1.22 

Beta 

.5059 

.335 

1.51 

90%  Confidence 

Limits 

Lower  limit 

Upper  limit 

Alpha 

-.7090 

.1012 

Beta 

-.0386 

1.0504 

R-squared  = 

.1872 

SC 

Standard  error  of 

estimate 

= 

.  8959 

GROUP=4 

AOV 

SOURCE 

DF 

SS  MS 

F 

Total 

11 

9. 

875 

Regression 

1 

1. 

848  1.848 

2.30 

Jan-Control 

1 

1. 

848  1.848 

2.30 

Residual 

10 

8. 

026  .803 

Variable 

Coefficient 

Standard  Error 

T 

Alpha 

-.2284 

.260 

-.88 

Beta 

.5447 

.  359 

1.52 

90% 

Confidence 

Limits 

Lower 

limit 

Upper  limit 

Alpha 

-.6618 

.2051 

Beta 

.5447 

1.1439 

R*squared  =  .2862 
Standard  error  of  estimate 


8949 


SC 

GROUP» 5 


ADV 


SOURCE 

DF 

SS 

MS 

F 

Total 

9 

8.976 

Regression 

1 

2.569 

2.569 

3.21 

Jan -Control 

1 

2.569 

2.569 

3.21 

Residual 

8 

6.407 

.801 

Variable 

Coefficient 

Standard  Error 

T 

Alpha 

-.1838 

.287 

-.64 

Beta 

.6573 

.  367 

1.79 

90%  Confidence  Limits 

Lower  limit 

Upper  limit 

Alpha 

-.6736 

.  3061 

Beta 

.0303 

1.2843 

TABLE  48 

ITERATIVELY  REWEIGHTED  LEAST 

SQUARES 

REGRESSION 

USING 

BI-WEIGHTS  FOR  THE  MODEL  OF  EQUATION 

GROUP  SIZE  3 

0=9 

A 

a 

A 

6 

Data  sets 

RN 

.048 

1.013 

FL 

-.063 

.664 

SC 

-.  116 

.  351 

C=4 

Data  sets 

RN 

.052 

1.031 

FL 

-.19  7 

1.049 

SC 

.154 

GROUP  SIZE  4 

.480 

0*9 

Data  sets 

RN 

-.  103 

.665 

FL 

-.093 

.  705 

SC 

-.169 

GROUP  SIZE  5 

.280 

0*9 

Data  sets 

RN 

-.065 

.865 

FL 

-.107 

.723 

SC 

-.126 

.413 
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2.  Maximum  Likelihood 


Table  49  displays  the  final  points  for  each  data  set 
along  with  the  inverse  Hessian  at  that  point.  Figures  103 
(RN) ,  104  (FL) ,  and  105  (SC)  are  interesting  in  that  they 
portray  the  contours  of  the  likelihood  functions  for  each 
data  set.  These  contours  show  the  variance  of  the  estimated 
parameters  in  a  graphic  way.  Note  how  data  sets  RN  and  FL 
seem  to  have  some  sort  of  horizontal  ridge  indicating  a 
good  pick  of  the  slope  parameter,  yet  the  contours  of  data 
set  SC  are  almost  circular  about  the  origin  of  the  axes 
indicating  no  significant  difference  from  zero  for  either 
parameter. 

TABLE  49 

MAXIMUM  LIKELIHOOD  ESTIMATES  OF  a  AND  3  ALONG  WITH 
ESTIMATES  OF  THEIR  VARIANCE  FOR  ALL  THREE  DATA  SETS 


A 


3 


a 


3 


a 

A 

3 


Data  set  RN 

.062 


2.918 


H 


-1 

L 


.257 

1.720 


Data  set  FL 
-.171 


.933 


.129 

.298 


Data  set  SC 
-.30  3 


.171 


.088 

.113 
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Figure  103.  Contours  of  log  likelihood 
function  for  data  set  RN 
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Figure  104.  Contours  of  log  likelihood 
function  for  data  set  FL 
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Figure  105.  Contours  of  log  likelihood 
function  for  data  set  SC 


C.  DISCUSSION 

A  recapitulation  of  all  the  parameter  values  is  in 
Table  50.  Some  interesting  observations  to  be  made  from 
this  table  are: 

(1)  The  slope  for  data  set  RN,  as  found  by  maximum 
likelihood,  is  much  greater  them  that  of  any  other  method 
or  data  set.  The  first  temptation  is  to  treat  this  as  an 
outlier,  yet  the  evidence  of  the  contour  plot  and  of  the 
validation  of  the  next  section  tend  to  back  up  this  number ■ 

The  reason  for  this  difference  is  a  possible  subject  for 
further  reserach. 

(2)  Except  for  the  intercept  values,  the  slopes  of  data 
sets  RN  and  FL  seem  to  be  fairly  consistent  within  and  between 
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regression  methods.  This  comment  is  made  in  light  of  the 
difference  of  these  data  sets  and  that  of  SC. 

(3)  Data  sets  RN  and  FL  seem  to  be  similar  in  many  ways, 
yet  data  set  SC  appears  to  be  different  in  both  degree  and 
signif icance . 

No  other  strong  pattern  is  apparent  in  these  parameter 
values.  Graphical  displays  of  the  parameters,  as  used  with 
grouped  data  are  given  in  Figures  106  (RN) ,  107  (FL) ,  and 
108  (SC).  Again,  note  the  significant  difference  of  the 
maximum  likelihood  line  for  data  set  RN.  Table  51  contains 
the  data  points  from  which  Figures  106,  107  and  108  were 
drawn. 

After  viewing  these  figures,  the  maximum  likelihood 
approach  is  the  preferred  method  for  the  Peninsula  data  sets, 
whereas  the  robust (C=4)  IRWLS  may  be  best  for  the  Valley  data 
set.  Although  fits  were  made  to  data  set  SC,  it  appears  as 
if  no  great  significance  has  been  found. 
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TABLE  50 

PARAMETER  FIT  RECAPITULATION 
FOR  ALL  DATA  SETS 


DATA  SETS 


RN 

FL 

SC 

Method 

1 .  Maximum 

a  : 

.062 

-.171 

-  .303 

Likelihood 

8  : 

2.918 

.933 

.171 

2.  OLS 

Group= 3 

a  : 

-  .195 

-.146 

-  .304 

8  : 

1.771 

1.052 

.506 

Groups 4 

a  : 

-  .170 

-.115 

-  .228 

8  : 

1.769 

1.357 

.545 

Group* 5 

a  : 

-  .265 

-.139 

-  .184 

8  : 

1.719 

1.523 

.657 

3.  IRWLS 

Groups  3 

a  : 

.048 

-.063 

-  .116 

C=9 

8  : 

1.013 

.664 

.351 

Groups 3 

a  j 

.052 

-.197 

.154 

C=4 

8  : 

1.031 

1.049 

.480 

Group  4 

a  : 

-.103 

-.093 

-  .169 

C=9 

8  : 

.665 

.705 

.280 

Groups 5 

a  : 

-.065 

-.107 

-  .126 

C=9 

8  : 

.855 

.723 

.413 
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PROBABILITY  OF  SUCCESS 


MODEL  FITS  TO  GROUP* 3 
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Figure  106.  Estimated  probability  of  greater-than- 
average  total  rest-of-year  rainfall  versus  the 
anomaly  of  logged  rainfall  for  January  for  data  set  RN 
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Figure  108 .  Estimated  probability  of  greater-than 
average  total  rest-of-year  rainfall  versus  the 
anomaly  of  logged  rainfall  for  January  for  data 
set  FL 
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TABLE  51a 


ACTUAL  VALUES  FOR  MODEL  FITS  OF  RN 


LOGGED 

ACTUAL 

MLE 

IRWLS 

IRWLS 

OLS 

ANOMALY 

VALUE 

C=9 

C=4 

-1.090 

0.000 

.042 

.258 

.255 

.107 

-  .460 

.  330 

.218 

.  397 

.396 

.267 

-  .190 

.  330 

.  379 

.464 

.464 

.  370 

.010 

.670 

.523 

.515 

.516 

.456 

.160 

.670 

.629 

.552 

.554 

.522 

.  360 

.670 

.753 

.602 

.604 

.609 

.500 

.670 

.821 

.635 

.638 

.666 

.980 

1.000 

.949 

.739 

.743 

.824 

TABLE 

51b 

ACTUAL  VALUES 

FOR  MODEL  FITS  OF 

FL 

LOGGED 

ACTUAL 

MLE 

IRWLS 

IRWLS 

OLS 

ANOMALY 

VALUE 

C=9 

C=*4 

-1.820 

.330 

.124 

.219 

.108 

.113 

-1.000 

0.000 

.249 

.326 

.223 

.232 

-  .460 

0.000 

.354 

.409 

.336 

.348 

-  .200 

.670 

.412 

.451 

.400 

.412  |i 

.110 

.  330 

.483 

.503 

.480 

.492  1 

.140 

.  330 

.490 

.507 

.487 

.500  1 

.200 

0.000 

.504 

.517 

.503 

.516  1 

.  340 

.670 

.5  36 

.541 

.540 

.553  } 

.594 

.500 

.670 

.573 

.567 

.581 

.610 

.670 

.598 

.585 

.609 

.621  J 

.750 

1.000 

.629 

.607 

.643 

.655  i 

1.140 

1.000 

.709 

.667 

.731 

.741  i 

TABLE  51c 


ACTUAL 

VALUES  FOR 

MODEL  FITS 

OF  SC 

LOGGED 

ACTUAL 

MEL 

IRWLS 

IRWLS 

ANOMALY 

VALUE 

C=9 

C=4 

-1.410 

.330 

.367 

.352 

.372 

-1.070 

0.000 

.  381 

.  380 

.411 

-  .960 

.330 

.385 

.389 

.424 

-  .550 

.330 

.402 

.423 

.473 

-  .260 

.670 

.414 

.448 

.507 

-  .100 

.330 

.421 

.462 

.526 

.020 

.670 

.426 

.473 

.541 

.200 

0.000 

.433 

.489 

.562 

.280 

.670 

.437 

.496 

.572 

.400 

.670 

.442 

.506 

.586 

.460 

1.000 

.444 

.511 

.593 

.560 

.670 

.448 

.520 

.604 

.690 

0.000 

.454 

.532 

.619 

.780 

.670 

.458 

.539 

.629 

.970 

.330 

.466 

.556 

.650 

1.250 

1.000 

.478 

.580 

.680 
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OLS 


.266 
.  300 
.312 
.  358 
.393 
.412 
.427 
.449 
.  460 
.475 
.482 
.495 
.511 
.523 
.547 
.581 


VII  .  VALIDATION  OF  LOGISTIC  MODELS 


A.  GENERAL 

The  various  parameters  that  were  estimated  in  the  previous 
section  all  may  be  subject  to  some  sort  of  validation.  How¬ 
ever,  this  paper  will  only  view  the  validation  for  the  maximum 
likelihood  approach  on  all  data  sets  and  the  IRWLS  (C*4) 
approach  on  data  set  SC.  The  validation  will  be  conducted 
against  the  reserved,  independent,  data  sets  of  years  1975 
through  1980.  These  are  the  same  data  sets  as  used  in  section 
IV. 

Table  52  portrays  the  reserved  data,  in  a  form  for  logistic 
analysis,  and  Figure  109  is  a  display  of  the  derived  contin¬ 
gency  tables  for  the  reserved  data. 
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TABLE  52 


RESERVED  DATA  IN  FORM  FOR  THE 
LOGISTIC  ANALYSIS 


DATA  SET  RN 


YEAR 

X 

COMPLEMENT 

Yh 

19  75 

-.739 

12.92 

o  fc 

1976 

-2.446 

11.28 

0 

1977 

-.477 

11.14 

0 

1978 

.883 

19.62 

1 

1979 

.542 

18.11 

1 

1980 

.752 

— 

— 

DATA  SET  FL 


1975 

-.470 

12.58 

0 

1976 

-2.610 

10.58 

0 

1977 

-.470 

11.10 

0 

1978 

.734 

19.32 

1 

1979 

.547 

13.72 

0 

1980 

.633 

— 

- 

DATA  SET  SC 


1975 

-.535 

17.97 

1 

1976 

-3.773 

11.19 

0 

1977 

.352 

10.03 

0 

1978 

1.054 

23.38 

1 

1979 

.512 

16.36 

0 

1980 

.405 

— 

- 
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Data  set  RN 


Control 


Complement 

-  3  1  0  |  3 

+  0  I  2  |  3 

3  I  2  |  5 


Control 


Data  set  FL 
Complement 

-  3  >  0  I  3 

+  1  I  1  I  2 

4  !  1  |  5 


Control 


Data  set  SC 
Complement 

-  1  1  *  I  2 

+  2  1  1  |  3 

3  j  2  I  5 


Figure  109.  2x2  contingency  Tables  of  reserved  data 

controlled  by  the  anomaly  of  January  rainfall.  The 
complement  is  the  anomaly  of  the  rest-of-year  rainfall. 
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B .  RESULTS 

1.  Data  set  RN 

The  model  proposed  by  the  maximum  likelihood 
parameters  is 


0t  " 


. 0618+2. 9183X. 
e  t 


VIII. 1 


l  +  e  t 

This  model,  when  applied  to  the  reserved  data  yields  Table  53. 

TABLE  53 

RESULTS  OF  LOGSTIC  VALIDATION 
ON  DATA  SET  RN 


YEAR 

X 

Y 

0 

1975 

-.739 

0 

.11 

1976 

-2. 746 

0 

o 

O 

• 

1977 

-.477 

0 

.21 

1978 

.883 

1 

.93 

19  79 

.542 

1 

.  84 

1980 

.752 

.91 

The  is  interpreted,  again,  as:  The  conditional 

probability  that  the  complement,  the  total  rainfall  for 
February  through  December,  will  be  above  its  mean  value, 
given  that  the  logged  January  anomaly  was  Xfc.  Thus  it 
appears  that  this  model  tends  to  predict  the  direction  of 
the  complements  deviation  well.  Figure  110  is  a  plot  of  the 
estimated  probabilities  against  the  actual  complement  anomaly. 
For  an  acceptable  fit,  this  plot  should  show  an  upward  to  the 
right  slope,  which  it  does. 
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Figure  110.  Plot  of  0^  versus  complement 
anomalies  for  data  set  RN 


2 .  Data  Set  FL 

This  data  set  is  quite  similar  to  the  RN  data, 
except  that  the  slope  parameter  is  only  a  third  of  that 
of  RN.  The  model  is 


171+ . 9325X. 

a  a  _£ _ t 

t  ,  7  7-7 171+. 9325X. 
l  +  e  t 


VI  I..2 
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TABLE  54 

RESULTS 

OF  LOGISTIC 

VALIDATION 

ON  DATA  SET 

FL 

Year 

X 

Y 

9 

1975 

-.470 

0 

.  35 

1976 

-2.610 

0 

.07 

1977 

-.470 

0 

.35 

1978 

.734 

1 

.63 

1979 

.547 

0 

.58 

1980 

.633 

- 

.60 

A  plot  of  the  probabilities  against  the  complement  anomalies 
is  in  Figure  111. 
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Figure  111.  Plot  of  9fc  versus  complement 
anomalies  for  data  set  FL 

This  fit  is  not  as  good  as  that  for  data  set  RN. 

The  outlier,  or  false  prediction  of  1979  may  not,  however, 
be  far  out  of  line.  The  sparsity  of  points  for  which  the 
complement  anomaly  was  positive  detracts  from  the  validation 


effort. 
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VII.  3 


VII.  4 


TABLE  55 


RESULTS  OF  VALIDATION  ON  DATA  SET  SC 


Year 

X 

Y 

e(D 

e(2> 

1975 

-.535 

1 

.40 

.47 

1976 

-3.773 

0 

.28 

.16 

1977 

.352 

0 

.44 

.58 

1978 

1.054 

1 

.47 

.66 

1979 

.512 

0 

.45 

.60 

1980 

.405 

- 

.44 

.59 

and  the  plot  of  the  probabilities  versus  the  complement 
anomalies  is  in  Figure  112. 
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Figure  112.  Plot  of  9fc  versus  complement 

anomalies  for  maximum  likelihood  and  IRWLS 
parameters  for  data  set  SC 

C.  DISCUSSION 

The  validation  of  the  maximum  likelihood  models  for  data 
sets  RN  and  FL  appears  to  be  acceptable.  However,  data  set 
SC  does  not  appear  to  be  acceptably  modeled.  In  fact,  as 
Figure  112  shows,  the  complement  appears  to  be  almost  inde¬ 
pendent  of  the  control.  This  is  also  shown  by  Figure  105 
where  it  can  be  seen  that  the  contours  are  very  flat  and 
circular  about  the  origin  of  the  (a, 8)  coordinate  system. 

The  one  apparent  outlier  of  data  set  FL  may  be  viewed 
as  very  close,  therefore  that  model  can  also  be  assumed 
to  be  validated. 


i 


i 
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VIII.  FURTHER  FINDINGS 


A.  SUMMER  MONTHS 

The  further  investigation  of  the  summer  months,  to 
parallel  the  modeling  of  the  winter  months,  yielded  some 
interesting  results.  These  results  are  shown  here  with 
no  attempt  at  analysis. 

The  summer  months  appear  to  be  increasing  in  total 
rainfall  and  in  variance.  This  is  more  true  for  the 
Peninsula  data  sets  than  for  the  Valley  data  set.  Figures 
113  (RN) ,  114  (FL) ,  and  115  (SC)  show  the  by-month  series 
of  summer  months.  The  total  summer  rainfall  series  by 
year  are  shown  in  Figures  116  (RN) ,  117  (FL) ,  and  118  (SC) . 
The  reserved  data  are  not  included,  yet  it  can  be  shown  to 
continue  the  indicated  trends. 
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RAINFALL  8N0MHL  ILS  IN  INCHES,  ‘  ^  _  ki  RAINFALL  ANOMALIES  IN  INCHES 


DATA  IS  FROM  'IE..PN 


:  113. 
means 


Monthly  plot  of  summer  months 
removed,  for  data  set  RN 


Figure  114.  Monthly  plot  of  summer  months 
only,  means  removed,  for  data  set  FL 
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DAT*  15  FROM  P, Erl 50 


SUMMER  MONTHS  OF  DRTfl 


Figure  115.  Monthly  plot  of  summer  months 
only,  means  removed,  for  data  set  SC 


DflTri  13  FROM  NCUPN 


Figure  116.  Yearly  plot  of  total  summer 
rainfall  for  data  set  RN 
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SUMMERS  TOTAL  3  r  rESR 


Figure  117.  Yearly  plot  of  summer  month 
rainfall  for  data  set  FL 
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Figure  118.  Yearly  plot  of  total  summer 
rainfall  for  data  set  SC 


A.  SIGNIFICANCE  OF  JANUARY 


The  identification  of  January  as  a  possible  predictor 
for  its  eleven  month  complement  raises  further  questions. 
One  of  the  questions  is  in  determining  which  part  of  the 
complement  lends  the  most  towards  its  predictability. 

Figure  119  is  a  plot  of  the  log-odds  and  chi-square  statis¬ 
tics  for  the  cumulative  complements.  Each  progressive 
column,  to  the  right,  of  Figure  119  indicates  these  statis¬ 
tics  for  another  cumulated  month,  i.e.,  the  first  column 
compares  the  anomalies  of  January  and  February  by  itself, 
the  second  coluim  is  a  comparison  of  January  to  February 
plus  March,  and  so  on  until  the  last  column  is  a  comparison 
of  January  to  the  entire  eleven  month  complement. 

Several  occurrences  to  be  noted  from  the  figure  are: 

(1)  The  log-odds  are  consistently  greater  than  zero. 

(2)  The  lack  of  increased  odds  and  significance  during 
the  summer  months. 

(3)  The  similarity  of  RN  to  FL  and  their  combined 
difference  to  SC  in  the  fall. 

These  indications  suggested  a  further  look  at  January 
versus  the  fall  months  only.  This  analysis  is  displayed 
in  Figure  120.  The  vertical  scales  of  Figure  119  and 
Figure  120  are  the  same,  yet  the  horizontal  scales  differ. 
This  figure  has  five  major  divisions.  The  left-most 
division  looks  at  January  versus  singular  months  in  the 
fall.  The  second  division  looks  at  January  versus  pairs 
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of  months  in  the  fall,  and  so  on  until  the  right -most 
column,  which  is  January  versus  the  total  fall  rainfall. 

This  figure  yields  no  apparent  significance,  and  unstable 
odds. 

The  combined  information  of  Figures  119  and  120  are 
mildly  confusing.  One  possible  explanation  may  be  that 
the  summer  months  somehow  cumulate  significance  and  devia¬ 
tion  direction,  in  order  to  allow  the  fall  contribution. 

This  possible  synergistic  affect  should  be  explored  further. 
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LOG-ODDS 


IX.  SUMMARY 


The  analysis  of  rainfall  data  is  carried  out  in  a 
comprehensive  way.  The  autoregressive  Markovian  model 
of  the  early  sections  could  not  stand  up  to  validation, 
but  it  did  point  to  some  sort  of  dichotomy  between  the 
seasons . 

2x2  contingency  analysis  was  effective  in  that  it 
brought  attention  to  the  predictive  ability  of  January. 
This  identification  of  January,  when  followed  by  the 
logistic  analysis  was  seen  to  be  successful  in  two  of 
the  three  data  sets.  Thus,  the  primary  conclusion  of  this 
thesis  is  the  predictive  ability  of  January  rainfall. 

The  physical  reasoning  behind  this  finding  must  be 
left  to  the  meteorologist.  Further  study  of  the  approach 
used  here  may  lead  to  improvement  in  seasonal  or  annual 
rainfall  forecasts  for  certain  climatic  regions. 
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APPENDIX  A 


DATA  SET  RN 


RAW 


YEAR 

OCT 

APR 

NOV 

MAY 

SEC 

JUN 

JAN 

JUC 

FEB 

AUG 

MAR 

SEP 

1931 

.  900 

1 .  100 

2.  780 
.  160 

6.330 

.240 

10.040 

.060 

2.980 

.080 

4.410 

.030 

1932 

.  170 

1 . 670 

2. 190 
.  490 

5.900 

.220 

2.060 

O.OO0 

.040 
.  110 

1 . 200 
.070 

1953 

.  380 
.350 

2.  130 
.420 

.360 

.400 

4.260 

0.000 

2.268 
.  160 

4.  910 
.030 

1934 

.  030 

1 . 740 

2.500 
.  370 

3.  130 
.060 

5.820 

0.008 

t .  740 
0.000 

.  160 
0. 000 

1953 

.830 

1.670 

1.960 

.390 

9.790 

0.000 

6.090 
.  070 

.  068 
.028 

.  130 
.360 

1936 

1 . 000 
1.490 

0.000 
2.  390 

.  840 
.200 

4.650 

0. 000 

3.328 

.070 

1.920 
.  030 

1937 

1.380 

4.710 

.930 

.360 

3.700 

.330 

3.710 
.  840 

3.660 

8.000 

7.  170 
.480 

1933 

.040 

.290 

.310 
.  120 

.490 

0.000 

4.330 

0. 000 

3.760 
.  040 

.320 
3.  140 

1939 

0.000 

.890 

0.000 

.340 

.390 

0.000 

4.300 
.  038 

4.338 

3.  000 

.  340 
.  130 

1960 

.  070 
1.290 

2.060 
.  720 

.950 

0.000 

1.890 

0.000 

1 .  170 
.  140 

2.380 

.090 

1961 

.  040 
.300 

1 . 740 
.  130 

1.190 

.230 

2.690 

0.000 

5.  170 
.230 

2.578 
.  130 

1962 

1 . 330 
3.930 

.373 

.660 

2.210 

.040 

3.050 

.040 

2.  700 
.001 

4.  140 
.400 

1963 

1.460 

.220 

3.770 

.060 

.330 

.220 

3.300 

.090 

.  420 
.338 

2.230 

.010 

1964 

.730 

2.260 

3.290 

.170 

6.430 
.  130 

2.560 

.030 

1.050 
.  160 

2.  440 
.020 

1965 

.230 

.270 

6.490 
.  130 

3.  560 
.120 

2.  320 
.290 

1.800 

.090 

.430 

.320 

1966 

.090 

7.  110 

4. 740 

.400 

4.  180 
1.360 

5.290 

.020 

.450 

.060 

3.480 
.  170 
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YEAR 

OCT 

APR 

NOV 

NAY 

1967 

.300 

.790 

1.610 

.320 

I960 

.310 

2.700 

3.  130 
.  120 

1969 

.500 

.330 

.720 

.050 

1970 

.590 

1.190 

6.  170 
.710 

1971 

.090 

.880 

1.990 

.090 

1972 

2.460 
.  130 

3.930 

.060 

1973 

2.208 

3.400 

3.870 

.030 

1974 

1.540 

1.760 

.560 

.010 

DEC 

JUN 

JAN 

JUL 

FEB 

AUG 

2. 270 
.010 

3.  100 
.060 

1.400 

.230 

3.270 

.  420 

9.450 

.040 

7.310 

.001 

3.080 

.300 

5. 910 
.030 

2.040 

.060 

4.990 

.030 

1.080 

.070 

.620 
.  130 

4.760 
.  150 

1.230 

.060 

1.030 

.040 

2.080 

.020 

6.  050 
.020 

3.880 

.030 

4.  730 
.  3  70 

3.  730 
.250 

.910 

.020 

2.488 

.170 

1.340 
.  170 

3.620- 

.430 

MAR 

SEP 

s.e-se 

.050 

1.310 
.  120 

2.  970 

.020 

1 . 360 
.430 

.030 
.  100 

4.520 
.  340 

4.460 

.010 

4.060 

.020 
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Logged  Anomalies 


YEAR 

OCT 

RPR 

NOV 

MAY 

SEC 

JUN 

JAN 

JUL 

FEB 

AUG 

MAR 

SEP 

1931 

.204 

112 

.  237 
-.  170 

.753 

.040 

.  363 
.003 

.272 

-.017 

.556 
-.  108 

1932 

-.231 
.  128 

.063 

.080 

.  620 
.023 

-.421 

-.034 

-1 . 065 
.010 

-.  344 
-.118 

1953 

-.116 

-.239 

.033 
.  032 

-.367 
.  161 

.  121 
-.054 

.  078 
.034 

.  644 
-.  137 

1934 

-.409 
.  134 

.  160 
.  307 

.  106 
-.117 

.331 

-.034 

-.096 

-.094 

- .  984 
-.185 

1935 

-.361 
.  128 

-.007 
.  143 

1 . 067 
-.  176 

.420 

.014 

-1 . 046 
-.074 

-.993 
.  122 

1 936 

.233 

.058 

-1.092 

.902 

-.702 

.007 

.  193 
-.054 

.404 
— .  026 

-.061 
-.  136 

1937 

.310 

.888 

-.435 
.  126 

.236 
.  125 

.011 

-.014 

.  792 
-.094 

.968 

.207 

1953 

-.399 

-.399 

-.680 

-.203 

-.913 

-.176 

.228 

-.034 

.807 

-.035 

-.833 

1.235 

1959 

-.438 

-.223 

-1.092 

-.026 

-.348 

-.176 

.  129 
-.024 

.606 

-.094 

-.323 

-.063 

1960 

-.370 

-.025 

.026 

.224 

-.697 
-.  176 

-.478 

-.034 

-.329 

.037 

.  143 
-.099 

1961 

-.399 

-.392 

-.084 
-.  179 

-.528 

.031 

-.233 

-.054 

.715 
.  129 

.  140 
-.046 

1962 

.408 

.741 

-.778 
.  188 

-.  146 
-.  136 

-.  140 
-.014 

.204 

-.893 

.504 
.  131 

1963 

.462 

-.633 

.470 

.302 

-.387 

.023 

-.035 

.032 

-.754 

.206 

.040 

-.175 

1964 

.  139 
.328 

.364 
-.  162 

.696 

-.036 

-.269 

-.005 

-.336 

.034 

.  103 
-.  166 

1963 

-.231 

-.613 

.921 
-.  196 

.369 

-.062 

-.339 
.  193 

-.046 

-.008 

-.773 

.092 

1966 

-.332 

1.239 

.633 

.018 

.333 
.  764 

.380 

-.034 

-.733 

-.036 

.  736 
-.023 
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YEAR 

OCT 

RPR 

NOV 

MAY 

DEC 

JUN 

JRN 

JUL 

FEB 

RUG 

MAR 

SEP 

196? 

-.116 

-.272 

-.  133 
-.041 

-.127 

-.166 

-.  123 
.005 

-.229 

.113 

.269 
-.  13? 

1969 

-.  169 
.454 

.326 

-.205 

.  140 
.175 

.  808 
-.014 

1.013 

-.093 

-.295 

-.072 

1969 

-.033 

-.554 

-.550 

-.270 

.  094 
.  087 

.  394 
-.024 

.008 

-.036 

.246 
-.  166 

1970 

.026 

-.070 

.878 

.218 

.478 
-.  146 

-.  80? 
.014 

-.622 
.  023 

-.047 
.  172 

1971 

-.352 

-.223 

.003 

-.232 

,439 

-.036 

-.737 

.005 

-.336 

-.055 

-1 .  103 
-.090 

1972 

.903 

-.732 

.346 

-.260 

-.  187 
-.  156 

.414 

-.034 

.824 

-.045 

.576 
.  107 

1973 

.725 

.629 

.491 

-.289 

.434 
.  139 

.015 
.  169 

-.45? 

-.074 

.568 
-.  175 

1974 

.494 
.  161 

-.648 

-.309 

-.065 

-.019 

-.689 
.  103 

.426 

.264 

.489 

-.166 
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Reserved  Data 


YEAR 

OCT 

APR 

NOV 

MAY 

DEC 

JUN 

JPN 

JUt- 

FEB 

PUG 

MPR 

SEP 

1974 

1.543 

1 . 763 

.560 

.010 

2.430 
.  170 

1.34  0 
.  170 

3.620 

.430 

4.063 
.  028 

1975 

1 .  700 

1. 740 

.  520 
.  070 

.370 
.  170 

.  130 
.020 

2.  970 
.970 

1 . 520 
.  420 

1976 

.600 

.040 

.  720 
1.210 

2.030 
.  080 

1.740 
.  030 

.  330 
.  020 

1.750 
.  650 

1977 

.  140 
5.430 

.  540 
.020 

5.350 
.  080 

6.730 

.040 

4.  730 
.  001 

5.240 
.  290 

1973 

.020 

.580 

2.  130 
.290 

1 . 590 
.020 

4. 320 
.350 

4.520 
.  090 

4.410 

.020 

1979 

1 . 300 
1.770 

2.350 

.570 

3.  180 
.040 

5.950 
.  730 

4.730 
.  090 

2.400 
.  890 
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APPENDIX  B 


DATA  SET  FL 


RAW 


YEAR  OCT  NOV 

APR  MAY 


193? 

.640 

1.270 

2.250 

.  030 

1  930 

1 . 230 

.  770 

.  280 

.  630 

1939 

.  750 

.  480 

.  720 

.  130 

1940 

0.000 

.  790 

3.760 

.  470 

1941 

.300 

.280 

3.  660 

.  790 

1942 

1.130 

1.470 

1.300 

.  100 

1943 

.620 

.  420 

1.160 

.830 

1944 

1.130 

4.630 

.430 

.300 

1943 

1.370 

1.400 

.040 

.640 

1946 

.260 

4.310 

.330 

.210 

1947 

1 .  160 

.  470 

3.230 

.400 

1943 

2.  130 

.  330 

.080 

.420 

1949 

.030 

1 .690 

1.470 

.  190 

1930 

1. 760 

3. 000 

1.060 

.270 

1931 

.670 

3.330 

1.100 

.  120 

1932 

.  140 

2.210 

1.230 

.370 

DEC 

JUN 

JAN 

JUL 

FEB 

AUG 

MAR 

SEP 

4 . 600 
0.800 

4.990 
.  130 

7.398 

3. 000 

3.838 
.  170 

2.890 
.  070 

3.190 

0. 888 

1 . 628 

0. 300 

2.3?0 
.  138 

1.478 

8.808 

3.430 

0.000 

3. 280 
8.808 

4.248 
8. 888 

4.320 

.040 

6.  130 
8.000 

7.730 

0.080 

3. 308 
8 . 800 

7.378 

0.008 

3.430 
.  120 

2.910 
.  188 

2.430 
0. 000 

.970 

3.080 

3.  050 
8.000 

2.730 

.060 

4.240 
.  820 

3.338 

.240 

3.690 
.  100 

3.880 
.  100 

.710 
0. 000 

1 . 940 
.118 

.  930 
0.000 

3.130 

.068 

2.338 
.  830 

4.230 

3.000 

1.080 

0.000 

2.360 

0.000 

3.000 
.  130 

1.970 
.  120 

.  440 
0.000 

1.130 

8.000 

1 . 938 
0. 800 

1 . 940 
.  170 

.  090 
0.000 

2.140 

0.000 

4.510 

0.080 

3.200 
.  058 

l  .  700 
.080 

2.340 
.  170 

4. 730 
0.  000 

1.370 

.110 

3.  120 
8.000 

1 . 740 
.070 

1.640 
.  090 

2.370 
.  120 

2.068 

.110 

2.  190 
.  100 

’  .  400 
.  030 

6.040 
.  170 

3.810 

.040 

2.110 

.070 

4,350 
.  198 

3.040 
.  220 

1.710 

0.000 

0. 080 
.030 

.  778 
.  070 
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YEAR 

OCT 

APR 

NOV 

hay 

DEC 

JUN 

JAN 

JUL 

FEB 

AUG 

MAR 

SEP 

1953 

.280 
.  780 

2.05O 

.390 

.310 

.510 

3.330 

.020 

1.530 
.  160 

4.410 

.030 

1954 

.020 

1.660 

2.  170 
.960 

2.  820 
.050 

4.560 
.  060 

1.560 
.  100 

.  130 
.050 

1955 

.040 

1.560 

1.650 

.590 

8.310 

0.000 

4.320 

.060 

1.740 

.090 

.  180 
.220 

1956 

.650 

1.340 

0.000 

2.640 

.  770 
.210 

4.  320 
.  020 

2.610 

.050 

1.390 

.220 

1957 

2.020 

3.980 

.850 

.498 

3.380 

.220 

2.620 
.  158 

4.620 
.  100 

6.760 

.550 

1958 

.040 

.460 

.390 
.  120 

.400 

.020 

4.470 

.030 

5.630 

.090 

.340 
2.  170 

1959 

.030 

.870 

0.000 

.410 

.620 

.010 

3.990 

.080 

3.279 

.040 

.650 

.010 

1968 

.070 

1.200 

1.450 

.750 

.640 

.040 

1.840 

.040 

1.  110 
.  160 

2.420 
.  140 

1961 

.  100 
.270 

2.030 
.  150 

1.380 

0.000 

2.230 
.  130 

5.830 

.280 

2.400 

.080 

1962 

1.360 

3.320 

.  440 
.650 

2.268 

.050 

3.830 

.110 

2.270 

.040 

4.328 

.578 

1963 

1.360 
.  100 

5.260 

.698 

.540 

.190 

3.  190 
.  190 

.390 

.330 

2.380 

.050 

1964 

1.060 

2. 180 

3.  160 
.  140 

5.970 

.060 

2.370 

.060 

.930 
.  190 

2.430 

.010 

1965 

.  160 
.  170 

6.250 
.  140 

6.230 
.  190 

3.050 

.320 

1.490 

.200 

.510 

.410 

1966 

.  140 
6.810 

4.690 

.790 

3.800 

.890 

5.750 

.060 

.490 
.  160 

5.070 

.230 

1967 

.270 

.720 

1.010 

.290 

2.060 
.  040 

3.110 

0.000 

1.350 

.270 

3.  180 

.  100 

1968 

.430 

2.490 

3.300 
.  190 

2.750 

.260 

9. 480 
.  ISO 

7.960 

.240 

1.200 
.  190 
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i 

i 


YEAR 

OCT 

APR 

NOV 

MAY 

DEC 

JUN 

JAN 

JUL 

PE  8 

AUG 

MAR 

SEP 

1 969 

.410 

.430 

1.010 

.000 

2.800 

.320 

3.330 

.030 

3.020 
.  190 

1.360 

.040 

1970 

.340 

1.190 

7.440 

.300 

4.450 

.030 

.370 
.  110 

.  *>40 
.  170 

1.940 

.350 

1971 

.  170 
.  720 

1.350 
.  190 

4.680 

.290 

1.000 

.090 

.  440 
.  100 

.080 
.  100 

1972 

2.340 

.310 

3.320 
.  100 

1.830 

.030 

3.530 

.070 

6.070 

.220 

3.330 

.410 

1973 

1.900 

2.990 

4.730 
.  110 

3.270 

.810 

3.  788 
.410 

1.190 

.210 

4.300 
.  100 

>»*> 


Logged  Anomalies 


YEAR 

OCT 

RPR 

NOV 

MRY 

DEC 

JUN 

JRN 

JUL 

FEB 

RUG 

MRR 

SEP 

1937 

.911 

.383 

-.  187 
-.289 

.434 
-.  138 

.392 

.995 

.932 
-.  186 

.788 

.811 

1939 

.348 

-.349 

-.436 
.  131 

-.  141 
-.862 

.933 

-.071 

-.233 
-.  186 

.831 

-.024 

1939 

.976 

-.234 

-.613 

-.216 

-.363 
-.  139 

.  346 
-.871 

1.081 

-.186 

.473 

-.146 

1949 

-.484 

1.113 

-.423 

.947 

.439 

-.991 

.369 

-.971 

.948 

-.186 

.688 
-.  146 

1941 

.  194 
.743 

-.769 

.244 

.914 
-.  138 

.  182 
.942 

.  145 
.968 

.849 
-.  146 

1942 

.272 

.937 

-.  193 
-.243 

-.391 
-.  138 

.881 

-.871 

.  183 
-.847 

.473 
-.  126 

1943 

-.991 

-.926 

-.636 

.266 

.386 

.993 

.  148 
.924 

.988 

-.818 

-.647 
-.  146 

1944 

.272 

-.439 

.723 

-.876 

-.  191 
-.823 

-.713 

-.971 

.212 
-.  947 

.  139 
-.117 

1943 

.379 

-.737 

-.  132 
.  137 

.395 
-.  139 

-.666 

-.971 

.831 
-.  196 

.293 

-.896 

1946 

-.253 

-.311 

.663 
-.  147 

-.181 

-.917 

-1.833 

-.971 

-.462 

-.196 

-.  182 
-.  146 

1947 

.286 

.631 

-.622 

-.991 

-.  191 
.827 

-1.312 

-.871 

-.974 
-.  196 

.341 
-.  146 

1949 

.637 

-.719 

-.797 

.913 

.  166 
-.891 

-.485 

.996 

.  127 
.931 

.362 
-.  146 

1949 

-.434 

.199 

-.918 

-.164 

-.496 

-.923 

.819 

-.871 

-.219 

-.938 

-.213 

-.869 

1939 

.332 

-.873 

.379 

-.899 

-.834 

-.917 

-.279 

.934 

-.959 

-.919 

-.273 

-.897 

1931 

.929 

-.934 

.463 

-.223 

.692 

.927 

.899 

-.932 

-.984 

-.938 

.539 

.929 

1932 

-.333 

.913 

.  139 
.113 

.329 

.969 

-.481 

-.871 

-1.218 

-.837 

-.613 

-.878 

194 


YEAR 

OCT 

NOV 

DEC 

JAN 

FEB 

MAR 

RPR 

HAY 

JUN 

JUL 

AUC 

SEP 

1933 

-.237 

.  108 

-.999 

.863 

-.271 

.503 

-.219 

-.009 

.232 

-.831 

.343 

-.  11? 

1934 

-.464 

.  14? 

.071 

.318 

-.278 

-1.061 

.  182 

.335 

-.081 

-.013 

-.010 

-.09? 

1935 

-.444 

-.032 

.962 

.363 

-.210 

-1.018 

.  144 

.  126 

-.  130 

-.013 

-.019 

.033 

193  6 

.817 

-1.007 

-.698 

.274 

.063 

-.  122 

.834 

.934 

.061 

-.051 

-.057 

.053 

195? 

.622 

-.392 

.208 

-.111 

.308 

.865 

.818 

.061 

.869 

.069 

-.010 

.292 

1958 

-.444 

-.673 

-.933 

.301 

.673 

-.891 

-.417 

-.225 

-.110 

-.041 

-.019 

1.009 

1939 

-.434 

-1.807 

-.787 

.210 

.233 

-.683 

-.170 

.006 

-.  120 

.006 

-.066 

-.  136 

I960 

-.416 

-.111 

-.774 

-.334 

-.472 

.046 

-.807 

.222 

-.091 

-.032 

.043 

-.013 

1961 

-.388 

.  102 

-.402 

-.223 

.783 

.040 

-.337 

-.190 

-.  130 

.031 

.  141 

-.069 

1962 

.373 

-.642 

-.087 

.  177 

-.034 

.488 

.667 

.  163 

-.081 

.034 

-.066 

.303 

1963 

.373 

.827 

-.83? 

.035 

-.889 

.010 

-.701 

.  187 

.044 

.  103 

.180 

-.097 

1964 

.239. 

.418 

.673 

-.  183 

-.561 

.049 

.361 

-.287 

-.072 

-.013 

.068 

-.  136 

1963 

-.335 

.974 

.709 

.001 

-.306 

-.771 

-.639 

-.207 

.044 

.207 

.077 

.  190 

1966 

-.333 

.732 

.316 

.312 

-.820 

.620 

1.260 

.244 

.307 

-.013 

.043 

.061 

196? 

-.243 

.026 

-.  131 

.016 

-.364 

.247 

-.234 

-.003 

-.091 

-.071 

.  133 

.019 

1968 

-.  126 

.432 

.033 

.831 

.974 

-.393 

.434 

-.  164 

.  101 

.069 

.110 

.028 

195 


i 


YEAR 

OCT 

APR 

NOV 

MAY 

DEC 

JUN 

JAN 

JUL 

FEB 

AUG 

MAR 

SEP 

1969 

140 

-.438 

-.  309 
-.261 

.066 
.  148 

.482 

-.822 

.  173 
.068 

-.  133 
-.  107 

1970 

-.052 

-.012 

1.  126 
.  120 

.427 

-.081 

-.  772 
.034 

-.724 

.051 

-.  105 
.  154 

1971 

-.327 

-.254 

-.871 
-.  172 

.463 

.125 

-.705 

.015 

1  » 

G»  CO 

u 
o  * 

-1. 10? 
-.051 

1972 

.  730 
-.526 

.868 
-.  189 

-.222 
-.  100 

.482 

-.003 

.738 

.093 

.32? 

.  190 

1973 

.602 

.588 

.742 

-.234 

.  183 
.463 

.  162 
.273 

-.434 

.085 

.521 

-.051 

Reserved  Data 


YEAR 

OCT 

APR 

NOV 

NAY 

DEC 

JUN 

JAN 

Jut 

FEB 

AUG 

MAR 

SEP 

1974 

1.960 

1.380 

.720 
.  1  10 

2.430 
.  180 

1.700 

.380 

2.760 

.430 

4.680 

.080 

1973 

1.480 

1.660 

.470 

.070 

.410 
.  130 

.200 
.  140 

2.820 

1.330 

1.060 

.820 

1976 

.110 

.038 

.780 

.960 

1.660 

.110 

1.700 

.020 

.630 

.040 

1.900 

1.050 

1977 

.  170 
3.210 

.620 

.030 

3.370 
.  100 

3.780 

.090 

4.710 

.030 

3.  160 
.310 

1978 

.  iao 

.360 

2.220 
.  130 

1.340 

.010 

4.700 

.330 

3.630 
.  120 

3.710 

.040 

1979 

1.380 

1.520 

3.360 

.380 

3.830 

.048 

5. 120 
.900 

4.420 
.  140 

2.030 
.  130 

197 


APPENDIX  C 


DATA  SET  SC 


RAW 


YEAR 

OCT 

MOV 

DEC 

JAN 

FEB 

MAR 

APR 

MAY 

/UN 

JUL 

AUC 

SEP 

1926 

.429 

3.360 

1.640 

2.440 

7.380 

1.390 

1.979 

.280 

.  108 

0.000 

0.000 

.  130 

192? 

.800 

2.  100 

2.  930 

.880 

2.430 

3.380 

.960 

.050 

0.000 

0.000 

0.000 

0.000 

1928 

.029 

3.  120 

3.410 

1.260 

1.370 

2.890 

.  1.290 

0.000 

1.230 

0.000 

0.000 

0.000 

1929 

0.090 

0.000 

.330 

3.370 

3.260 

4.230 

1.330 

1.610 

0.000 

0.  000 

0.000 

.030 

1938 

.030 

1.480 

.030 

4.260 

1.560 

.920 

.490 

.660 

.400 

0.000 

0.000 

0.000 

1931 

.020 

1.730 

10.260 

4.320 

4.690 

.300 

.  140 

.360 

0.000 

0.800 

0.000 

0.000 

1932 

0.000 

.  180 

3.120 

6.920 

.920 

1.370 

.330 

.820 

.040 

0.000 

0.000 

0.000 

1933 

1.120 

0.000 

8.230 

3.  150 

3.290 

0.000 

.130 

.850 

.720 

0.000 

0.000 

.060 

1934 

.  130 

2.890 

2.690 

7.000 

.770 

4.880 

3.910 

0.000 

0.000 

0.000 

.630 

0.000 

1933 

.370 

.810 

2.010 

2.660 

9.970 

1.680 

2.460 

.480 

.  340 

.430 

0.000 

0.009 

1936 

.430 

0.000 

4.  160 

4.620 

6.770 

7.730 

.320 

0.000 

.130 

0.000 

0.000 

0.000 

t*3? 

.060 

1.010 

6.390 

3.440 

13.029 

3.099 

2.630 

0.000 

9.000 

0.000 

0.000 

0.000 

1933 

.810 

1.040 

2.480 

3.340 

2.70O 

3.000 

.420 

.220 

.  100 

0.000 

0.000 

.100 

1939 

1.040 

.3  70 

2.  380 

9.200 

0.960 

2.320 

.600 

.270 

0.000 

0.000 

0.000 

.310 

1948 

1.040 

.370 

2.380 

9.200 

8.960 

2.320 

.600 

.270 

0.000 

0.000 

0.000 

.310 

1941 

.430 

.330 

9.430 

9.  190 

10.420 

0.920 

3.000 

.610 

.  100 

0.000 

0.000 

0.000 

198 


YEAR 

OCT 

APR 

NOV 

MAY 

DEC 

;un 

JAN 

JUL 

PE8 

AUG 

1942 

1.130 
4.  910 

1.040 

1.010 

10.490 

0.000 

5.  710 
0.000 

2.460 
0. 009 

1943 

.740 

1.070 

1.960 

9.000 

1.630 

.090 

3.400 

0.000 

2.429 
0.  000 

1944 

.740 

1.940 

.360 

.758 

3.330 

.320 

3.940 

0.000 

9.239 

0.000 

1945 

1.030 

.280 

3.240 

.200 

2.390 

0.000 

1.330 

0.000 

7.610 

.240 

1946 

2.000 

0.000 

2.050 

.430 

7.690 

0.000 

.350 

0.000 

2.770 

0.000 

1947 

.790 

3.860 

.670 

.620 

2.360 
.  150 

.050 

0.009 

2.480 

0.000 

1948 

1.740 
.  150 

.  ISO 
.510 

6.010 

0.000 

1.180 

.020 

3.040 

.039 

1949 

.  190 
1.470 

1.470 

.260 

1.690 

0.000 

6.09O 

0.000 

2.970 

.050 

1950 

2.720 

1.370 

6.330 

.690 

3.040 

.050 

2.950 

0.098 

1.940 

0.000 

1951 

1.440 

.800 

4.010 

.220 

8.360 

0.000 

9.388 

8.808 

1.619 

0.000 

1952 

.070 

1.090 

3.000 

.620 

8.760 

.  100 

2.638 

8.888 

0.000 

.090 

1953 

.300- 

.530 

2.  160 
.  190 

.490 

.250 

4.378 

3.880 

3.280 

0.000 

1954 

0.000 

2.840 

1.980 

1.360 

2.980 

0.000 

5.530 

0.000 

1.990 

0.000 

1955 

.030 

1.600 

2.240 

.480 

14.940 

0.000 

6.590 

0.000 

2.250 

0.000 

1956 

.790 

1.720 

.020 

1.930 

.530 

.090 

5.210 

0.000 

5.110 

0.000 

1957 

1.470 

8.500 

.930 

.560 

3.990 

.140 

4.630 

0.000 

9.680 

0.000 

t 


MAR 

SEP 

3.280 

0.000 

3.740 

0.000 

1.000 

0.000 

4.380 

0.000 

3.280 

0.000 

4.380 
6. 000 

5.960 

0.000 

2.230 

0.000 

1.040 

.030 

6.820 

.070 

2.090 

0.000 

4.680 

0.000 

.310 

0.000 

.510 

.210 

1.570 

.070 

r.820 

.780 


199 


VE8* 

OCT 

8PR 

NOV 

MflV 

DEC 

JUN 

JflN 

JUL 

PE8 

flue 

tlftR 

SEP 

1938 

.038 

.  tie 

1.810 

.099 

.278 

8.808 

6.330 
9. 000 

6.920 

.090 

.238 
2.  100 

1959 

9.888 

1.918 

9.909 

.338 

.710 

0.808 

3.260 

8.000 

6.320 

8.908 

.740 

.038 

1988 

.858 

.868 

3.368 

.328 

1.810 

.078 

1.918 

9.088 

.910 

8.000 

2.490 

0.000 

1961 

.178 

.118 

2.628 

.349 

1.308 
9.  088 

2.908 

8.008 

11.290 

8.000 

2.960 

0.808 

1962 

1.318 

3.998 

.  108 
.238 

2.290 

.830 

7.818 

0.988 

2.208 

0.008 

4.420 

.890 

1963 

1.388 

.398 

3.648 

1.338 

.460 

.178 

4.988 

8.988 

.468 

.380 

3.888 

0.880 

1964 

1.298 

2.838 

3.399 

.148 

4.339 

.028 

2.838 

0.880 

1.038 

.230 

2.320 

8.080 

t963 

.878 

.278 

6.639 

9.888 

4.798 

.880 

1.318 

.290 

1.340 

8.008 

.688 
,  .180 

1966 

8.888 

6.418 

3.748 

.328 

3.418 

.190 

6.898 

0.088 

.438 

0.088 

6.  168 
.138 

196? 

.278 

.678 

1.738 

.438 

1.970 
.  128 

3.368 

9.088 

1.818 

.028 

3.298 

8.088 

1968 

.338 

2.  140 

1.938 
.  128 

4.068 
.  188 

13. 160 
8.088 

11.970 

0.990 

1.028 

8.888 

1969 

.340 

.010 

.798 

.078 

3.  110 

.048 

7.088 

0.988 

4.610 

8.808 

1.310 

0.088 

1978 

.  110 
1.210 

3.868 
.  138 

6.440 

8.000 

1.208 

0.088 

.538 

.848 

1.490 

.170 

1971 

.240 

.090 

1.438 

.888 

3.368 

.828 

1.198 

0.088 

.988 

9.888 

.010 

.030 

1972 

2.330 

.  100 

3.790 

.828 

2.668 

8.080 

7.390 

8.008 

9.288 

8.808 

4.430 

.070 

1973 

1.810 

1.970 

6.318 

8.888 

2.360 

.110 

4.440 

.878 

1.980 

0.808 

4.740 

e.eoe 

200 


Logged  Anomalies 


YEAR  OCT 
APR 

NOV 

NAY 

DEC 

JUN 

1926 

'.093 

-.090 

1.273 

-.073 

-.437 

.003 

192? 

.  144 
144 

.170 

-.271 

-.063 

-.092 

192  a 

-.424 

.811 

.433 

-.320 

.076 

.710 

1929 

-.444 

.029 

-.961 

.639 

-.304 

-.092 

1936 

-.393 

-.418 

-.033 
.  187 

-1.339 

.244 

1931 

-.424 

-.686 

.030 

-.013 

1.013 

-.092 

1932 

-.444 

-.332 

-.796 

.279 

.008 

-.033 

1933 

.388 

-.693 

-.961 

.293 

.81? 

.430 

1934 

-.384 

1.116 

.397 

-.320 

-.  102 
-.092 

1933 

-.  129 

.424 

-.368 

.072 

-.306 

.200 

1936 

-.086 

-.398 

-.961 
-.320  . 

.233 

.030 

193? 

-.306 

.472 

-.263 

-.320 

.619 

-.092 

1938 

.138 

-.466 

-.240 

-.121 

-.161 

.003 

1939 

.269 

-.34? 

-.646 

-.081 

-.  190 
-.892 

1948 

.269 

-.34? 

-.646 

-.081 

-.  190 
-.092 

1941 

-.072 

.966 

-.676 

.136 

.930 

.603 

JAN 

JUL 

FES 

AUG 

MAR 

SEP 

-.  34? 
-.013 

.785 

-.831 

-.372 

.063 

-.932 

-.013 

-.212 

-.031 

.279 

-.073 

-.768 

-.013 

-.381 

-.031 

.113 

-.073 

.269 

-.013 

.003 

-.031 

.413 

-.046 

.077 

-.013 

-.304 

-.031 

-.591 

-.073 

.089 

-.013 

.294 

-.031 

-.633 

-.073 

.406 

-.813 

-.792 

-.031 

-.299 

-.073 

-.160 

-.013 

.393 

-.031 

-1.243 

-.01? 

.496 

-.013 

-.873 

.469 

.329 

-.073 

-.283 

.343 

.931 

-.831 

-.237 

-.073 

.  143 
-.013 

.606 

-.031 

.929 

-.073 

-.092 

-.013 

1.196 

-.031 

.964 

-.073 

-.978 

-.013 

-.  136 
-.031 

.  163 
.020 

.739 

-.013 

.854 

-.031 

-.043 
.  195 

.739 

-.013 

.854 

-.031 

-.043 

.195 

.238 

-.013 

.991 

-.031 

1.052 

-.073 

201 


YEAR  OCT 
RPR 


1942 

.312 

.960 

1943 

.  110 
-.098 

1944 

.110 

.261 

1945 

.289 

-.370 

1946 

.891 

-.817 

1947 

.  136 
.764 

1948 

.364 

-.677 

1949 

-.270 

.087 

1930 

.870 

.046 

1931 

.448 

-.229 

1932 

-.376 

.244 

1933 

-.  191 
-.392 

1934 

-.444 

.328 

1933 

-.414 
.  138 

1936 

.  138 
.  183 

1937 

.460 

1.474 

NOV 

MAY 

DEC 

JUN 

-.248 

.378 

1.033 

-.092 

.  124 
-.320 

-.441 

-.006 

-.634 

.239 

.037 
.  183 

.483 
-.  138 

-.187 

-.092 

.  134 
.837 

.734 

-.892 

-.448 
.  162 

-.  196 

.048 

-.821 

.092 

.339 

-.892 

-.037 

-.009 

-.419 

-.092 

1.031 

.204 

-.012 

-.043 

.630 

-.121 

.820 

-.092 

.423 
.  162 

.370 

.003 

.  189 
-.  146 

-1.009 
.  131 

.131 

.538 

-.027 

-.092 

.214 

.872 

1.361 

-.092 

-.941 

.733 

-.903 

-.006 

.384 
.  124 

.  199 
.039 

JAN 

JUL 

FEB 

AUG 

.  321 
-.013 

-.20  3 
-.031 

.638 

-.015 

-.213 

-.031 

.024 

-.015 

.386 

-.031 

-.737 

-.013 

.  709 
.  184 

-.968 

-.013 

-.  117 
-.031 

-1.334 

-.015 

-.  19? 
-.031 

-.804 

.003 

-.048 

-.002 

.376 

-.013 

-.066 

.017 

-.209 

-.013 

-.366 

-.031 

.797 

-.013 

-.485 

-.031 

-.280 

-.013 

-1.444 

.033 

.098 

-.013 

.010 

-.031 

.293 

-.013 

-.349 

-.031 

.444 

-.013 

-.266 

-.031 

.243 

-.013 

.366 

-.031 

.  143 
*.013 

.924 

-.031 

NAR 

SEP 

.2X1 

-.073 

.313 

-.073 

-.530 

-.073 

.529 

-.373 

.211 

-.073 

.440 

-.073 

.697 

-.073 

-.071 

-.073 

-.330 

-.046 

.814 

-.007 

-.113 

-.073 

.494 

-.073 

-.973 

-.075 

-.831 

.116 

.299 

.007 

.934 

.302 


202 


YEAR 

OCT 

APR 

NOV 

MAY 

DEC 

JUN 

JAN 

JUL 

FEB 

AUG 

195$ 

-.414 

-.713 

-.263 

-.234 

-1.169 

-.092 

.412 

-.015 

.623 
.  .033 

1959 

-.444 

.251 

-.961 

-.035 

-.372 

-.092 

.251 

-.813 

.373 

-.031 

1960 

-.  395 
-.  197 

.620 

.898 

-.375 

-.025 

-.315 

-.013 

-.  797 
-.031 

1961 

-.28? 

-.713 

.325 

-.023 

-.573 

-.092 

-.484 

-.013 

1.064 

-.031 

1962 

.589 
.  790 

-.866 

-.113 

-.217 

-.043 

.  593 
-.013 

-.281 

-.031 

1963 

.389 

-.488 

.574 

.534 

-1.030 

.065 

.  026 
-.015 

-1.066 

.231 

1964 

.385 

.531 

.543 
-.  189 

.335 

-.072 

-.240 

-.013 

-.726 
.  192 

1965 

-.376 

-.578 

1.071 

-.320 

.  348 
-.015 

-.330 

.240 

-.594 

-.031 

1966 

-.444 

1.186 

.595 

.098 

.430 

.082 

.376 

-.813 

-1.087 

-.031 

196? 

1  1 

cj  ro 

©  © 
*  <JI 

.030 

.03? 

-.354 

.021 

-.110 

-.013 

-.746 

-.012 

1968 

-.  159 
.327 

.  121 
-.207 

.213 

.003 

1.200 

-.013 

1.118 

-.031 

1969 

151 

-.224 

-.379 

-.233 

.003 

-.033 

.  496 
-.015 

.280 

-.031 

1970 

-.339 

-.024 

.965 

-.180 

.399 

-.092 

-.794 

-.013 

-1.019 

.008 

1971 

-.229 
-.  181 

-.073 

-.243 

.  473 
-.072 

-.  799 
-.013 

-.761 

-.031 

1972 

.904 

-.722 

.934 

-.300 

-.111 

-.092 

.368 

-.013 

.878 

-.831 

1973 

.509 

.271 

1.028 

-.320 

-.  196 
.012 

.  1 1  1 
.053 

-.332 

-.031 

MAR 

SEP 

-1.036 

1.056 

689 

-.026 

.00? 

-.  0?5 

.  133 
-.075 

.44? 

.011 

.  163 
- .  0?5 

.015 
-.  0?5 

- .  ?24 
.090 

.726 

.04? 

.213 
- .  8?5 

-.  540 
-.  0?5 

-.323 

-.075 

-.331 

.082 

1.233 

-.046 

.449 

'.087 

.504 

-.075 
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Reserved  Data 


YEAR 

OCT 

APR 

1974 

1.340 

1.660 

1973 

1.310 

1.460 

1976 

1.450 
.  140 

1977 

.020 

3.600 

1378 

0.000 

.380 

1979 

1.010 

2.330 

NOV 

MAY 

DEC 

JUN 

.760 

.820 

3.370 
.  040 

.  730 
.030 

.330 
.  180 

.600 

.970 

1.830 

.080 

.400 

.030 

3.680 

0.000 

2.720 

.260 

t.200 

8.000 

2.  140 
.420 

3.870 

0.000 

JPN 

3UL 

FEB 

«uc 

2.840 

.070 

3.490 
.  120 

.080 

.001 

2.  330 
1.030 

2.430 

0.000 

.  460 
0. 000 

9.990 

0.000 

3.400 
0.  000 

3.810 
.  120 

4.  660 
0.000 

3.228 

.930 

10.890 

0.000 
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hrr 

SEP 

s.200 

0.00a 

1.430 

.  730 

2.000 

.480 

7.070 

.360 

3.020 

0.  000 

2.390 

0.000 
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